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I. INTRODUCTION 

In recent years, the study of New Physics (NP) that lies close to the electroweak (EW) energy scale 
has been motivated primarily by the hierarchy problem. However, it is possible that the correct solution 
to this problem or the detailed nature of EW symmetry breaking remain to be proposed. Experimental 
input, as expected from the LHC, is crucial. Furthermore, hints for new physics (NP) may have already 
emerged from the Tevatron. In this paper we are motivated by recent experimental anomalies at the Tevatron 
and the strong discovery potential at LHC to explore collider signatures of new physics (NP) sectors that 
are flavor symmetric. They will be taken to be invariant under the global flavor symmetry group Gp = 
U(3) Ur x U(3) Dr x U(3) Ql , or its subgroup H F = U(2) Ur x U(2) Dr x U(2) Ql x U(l) 3 (where the 
quarks of the first two families are in doublets of the corresponding SU(2) factors). The group Gf is the 
global symmetry of the Standard Model (SM) in the limit where one can neglect the Yukawa interactions 

C Y = Yu u R H T ia 2 Ql - Y D d R Q L + h.c. , (1) 

where Yu and Yd are the up and down quark Yukawa matrices, respectively. 

The NP sectors will contain scalar or vector fields that have masses < TeV and 0(1) couplings to 
quarks. At the same time they will be consistent with flavor changing neutral current (FCNC) constraints 
precisely because of their flavor structure, as long as the breaking of the flavor symmetries is sufficiently 
small. In the SM the top and bottom Yukawa couplings break the flavor group Gf to its Hp subgroup. We 
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take this breaking into account in our analysis of NP sectors that are initially Gf symmetric. The existence 
of the Hp symmetry at low energies protects these theories against dangerously large FCNC's, e.g., in 
neutral meson mixing. This protection satisfies the naturalness criteria of Glashow and Weinberg [ 1 ] and is 
not the result of simply tuning parameters. Note that NP models that have an approximate Hp symmetry 
are sometimes referred to as models of next to minimal flavor violation [2J. The breaking of Hp in the SM 
is due to the other quark masses and the CKM mixing angles, and is thus small. The precise mechanism 
by which Hp is broken in the NP sector will not be important when we explore flavor diagonal collider 
signatures. However, the nature of Hp breaking will be relevant to our discussion of low energy FCNC's, 
see below. 

Scalar and vector fields with dimension four Gf invariant direct couplings to quarks are limited in their 
allowed charge assignments by flavor and the SM gauge symmetry. There are only 14 different nontrivial 
flavor representations allowed in each case. In the case of Hp symmetric models the possible NP fields are 
conveniently classified in terms of these representations, with the understanding that they need not come 
in complete Gf multiplets. A systematic exploration of new flavor symmetric sectors is therefore feasible, 
either in general, or with the aim of explaining a particular anomaly. 1 

New flavor symmetric sectors that are perturbatively coupled to quarks are particularly interesting to 
consider as candidate explanations for Tevatron anomalies. In the first part of this paper, we focus our 
attention on two > 3a anomalies: (i) the CDF measurement of the ti forward backward asymmetry, Ap B , 
for mi t > 450 GeV [6] is 3Aa away from the NLO SM prediction. (A recent D0 analysis [7] does 
not observe a significant m t i dependence in the "folded" detector level asymmetries, but it appears to be 
consistent with the CDF detector level measurements within errors.) The inclusive ti forward backward 
asymmetry, averaged over the CDF semileptonic and hadronic ||8| ti decay samples and the recent 
D0 measurement f71. is ~ 3a from the NLO SM prediction; (ii) the like sign dimuon asymmetry measured 
by D0 is 3.9a away from the SM expectation ||9l|T0l. Each of these anomalies, if confirmed, points to a 
relatively low scale of NP with a significant coupling to quarks. We identify flavor symmetric models that 
have the potential to explain them either individually or simultaneously, and study related constraints. In the 
case of Gf symmetric models, under the assumption that the NP only couples to quarks, some hierarchies 
among these couplings would be required in order to consistently explain the A FB anomaly, e.g., due to 
the absence of dijet or ti resonances at the Tevatron and LHC. Thus, breaking of the Gf symmetry to Hp 
would be necessary. Alternatively, one could consider Hp symmetric models where the more constrained 



1 Color symmetric fermion content that mixes with the SM fermion fields is not as constrained in its allowed representations. 
Initial studies of vector-like fermions have also been undertaken in Refs. |4||5]. Of the models studied only two were natural in 
the Glashow- Weinberg sense 1 5 1. 
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quark couplings would simply be absent. 

It is possible that the above anomalies could be due to statistical fluctuations or underestimates of the- 
oretical or experimental errors. Even if this turns out to be the case, the models we explore in this paper 
are interesting in their own right, as they have strong discovery prospects at LHC. Again, this is because 
their flavor symmetric structures allow for sub-TeV NP mass scales. In the second part of this paper, we 
address the phenomenology of flavor symmetric sectors more generally. The global flavor symmetries we 
consider could be accidental, or they could be a remnant of the underlying mechanism generating the SM 
flavor structure (such as non-Abelian horizontal symmetries). We do not concern ourselves with the UV 
origin of these symmetries, but instead focus on the collider and low energy phenomenology of the new sec- 
tors. This approach is inspired by effective field theories (EFT), where one generally constructs all possible 
interactions consistent with the symmetries of interest. The analysis of flavor diagonal collider constraints 
and signatures can then be kept quite general, i.e., independent of the way Hp is broken, as already men- 
tioned. For simplicity, in attempts to explain a measured deviation from the SM we will only consider the 
phenomenology of single Gf multiplets (or the corresponding Hp multiplets), effectively assuming that 
there is a significant mass gap with other possible representations. Moreover, we will only consider their 
couplings to quarks. Note that more generally, these fields could couple to additional states transforming 
under Gf or Hp, possibly providing them with additional decay channels. 

The determination of low energy flavor physics constraints on flavor symmetric models generally re- 
quires the breaking of Hp to be specified. When determining these constraints we assume the Minimal 
Flavor Violation (MFV) hypothesis lfTTl[T2l . i.e., that all breaking of Hp is due to the SM Yukawas. This 
enforces maximal consistency with FCNC constraints through a symmetry principle, and allows us to ex- 
plore how low the NP mass scale can be. In all the models we consider, the new states can have EW scale 
masses. In MFV models that lead to class-2 operators (those that involve right handed fields) in the lan- 
guage of 1 13 ], the breaking of Hp can actually be orders of magnitude larger then assumed in MFV, while 
still obeying the FCNC bounds. 

The paper is organized as follows. In Section [IT] we list all vector and scalar representations of the form 



we have motivated, and write down in detail the vector field Lagrangians for two examples. In Section III 



we systematically discuss the potential of models of this form to explain the A l p B anomaly, the D0 dimuon 



anomaly, and related phenomenology. In Section IV we explore existing bounds on these models from LEP, 
electroweak precision data (EWPD), FCNCs and dijet studies at the LHC and Tevatron. In Section [V| we 



discuss additional LHC phenomenology. Finally, in Section VI we give our conclusions. Many details have 
been relegated to the Appendices. In Appendix[A]we list the details of flavor symmetric vector Lagrangians, 
in Appendix [B] we gives the details of 2 — > 2 scattering calculations and phenomenology, and in Appendix 
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Case 


SU(3) C SU(2) L 


U(1) Y U(3) Ur 


x U(3) Dr x U(3) c 


j L Couples to 


Is,o 


1,8 


1 





(1,1,1) 


d~R 7 P d R 


H S ,o 


1,8 


1 





(1,1,1) 


ur 7 p ur 


IHs.o 


1,8 


1 





(1,1,1) 


Ql 7 p Ql 


iv S)0 


1,8 


3 





(1,1,1) 


Ql1» Ql 


V S , G 


1,8 


1 





(1,8,1) 


d~R l* 1 d R 


VI S , 


1,8 


1 





(8,1,1) 


Ur 7^ Ur 


VII S , 


1,8 


1 


-1 


(3,3,1) 


d~R 7 P Ur 


VIII S , 


1,8 


1 





(1,1,8) 


Ql 7 p Ql 


IX S , 


1,8 


3 





(1,1,8) 


Ql 7 p Ql 


^3,6 


3,6 


2 


-1/6 


(1,3,3) 


d R r Qi 


Xi 3,6 


3,6 


2 


5/6 


(3,1,3) 


UR^Ql 



TABLE I: The flavor and gauge representations for vector fields that can couple directly to quarks through Gf sym- 
metric dimension four interactions without the insertion of a Yukawa matrix. Q C L denotes the right handed conjugate 
representation of the left handed SM doublet. 

[C]we give a detailed discussion of constraints from meson mixing amplitudes. 

II. G F SYMMETRIC REPRESENTATIONS 

We are interested in scalars and vectors that couple directly to quarks through dimension four interac- 
tions. The scalar fields of this form are renormalizable models, while the vector fields are nonrenormaliz- 
able. In this section we list all the possible representations of Gf and the SM gauge group that such fields 
can have when Gf is unbroken. The vector field representations are listed in Table [I] and the scalar field 
representations are listed in Table [TT] (the latter have been studied and classified in (51 Q31). This completes 
the program initiated in [3]. The complete set of Hp symmetric representations which can couple directly 
to quarks through dimension four interactions appear in these tables as submultiplets of the Gf representa- 
tions. For example, in model VI S ( Q ), the corresponding Hp symmetric vector representations would consist 
of a triplet, a complex doublet and a singlet of SU(2)u R , which are color singlets (octets), or a subset of 
these. 

Several remarks are in order before we construct the Lagrangians. 

• I Sj0 , II Sj0 and III s ,o carry the same quantum numbers and are thus a sub-classification of interac- 
tions of a single field. For instance, in the case of a color singlet vector with the same couplings to 
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Case 


SU(3) C SU(2) L 


U(1) Y U(3) Ur 


x U(3) Dr x U(3) Ql 


Couples to 


Si 


1 


2 


1/2 


(3,1,3) 


ur 


Ql 


Sii 


8 


2 


1/2 


(3,1,3) 


Ur 


Ql 


Sin 


1 


2 


-1/2 


(1,3,3) 


d R 


Ql 


Siv 


8 


2 


-1/2 


(1,3,3) 


d R 


Ql 


S v 


3 


1 


-4/3 


(3,1,1) 


Ur 


ur 


Svi 


6 


1 


-4/3 


(6,1,1) 


Ur 


Ur 


Svn 


3 


1 


2/3 


(1,3,1) 


d R 


d R 


Sviii 


6 


1 


2/3 


(1,6,1) 


d R 


d R 


Six 


3 


1 


-1/3 


(3,3,1) 


dn 


ur 


Sx 


6 


1 


-1/3 


(3,3,1) 


d R 


Ur 


Sxi 


3 




-1/3 


(1,1,6) 


Ql 


Ql 


Sxn 


6 




-1/3 


(1,1,3) 


Ql 


Ql 


Sxiii 


3 


3 


-1/3 


(1,1,3) 


Ql 


Ql 


Sxiv 


6 


3 


-1/3 


(1,1,6) 


Ql 


Ql 


Sh,8 


1,8 


2 


1/2 


(1,1,1) 


Qlur, 


QlAr 



TABLE II: Different scalar representations that are not singlets under the flavor group that are Gf symmetric [3 1 (the 
upper rows). The two flavor singlet representations are in the last row and were discussed in lfT4l . 

QL,UR,dji this is just the baryonic Z'. We found it useful to split the interactions into three sub- 
groups. At colliders there is no interference among these interactions up to effects suppressed by 
light quark masses (but if there are relations between their couplings this can have important con- 
sequences for the predicted cross section; for example, for a purely axial gluon the NP interference 
with the SM amplitude does not contribute to the top pair production cross section JT514T71 ). In the 
treatment of FCNCs the interference effects are trivial to include in the analysis. 

• Many of the scalar and vector fields do not lead to proton decay at any order in perturbation theory 
due to the SM gauge symmetry and Gf- The vectors X-XI and scalars Sy-Sxiv carry baryon number 
±2/3 and may lead to proton decay if they also couple to lepto-quark bilinears, e.g., L c l ^^ur and 
L^7 M dij for fields X and XI, respectively. This type of coupling is not possible for scalars or vectors 
in the color 6 representation and can be forbidden for the 3 color representation fields by extending 
the flavor group to the lepton sector of the SM . 

• We assume that the new quanta have weak scale masses and that the cut-off of the theory is well 
above the weak scale so that we only need to focus on dimension four interactions for most of our 
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discussion. Other dimension four couplings such as BpTV^P), with B^ u the hypercharge field 
strength and the vector fields, are not directly relevant to the phenomenology of interest in this 
paper. We leave the exploration of these interactions to a future publication. 

• The kinetic terms (with flavor breaking insertions) can always be made universal through field redef- 
initions. Below, we only write down the interaction terms. 

• Tree-level exchanges of fields in a single representation of Gf cannot explain both of the Tevatron 
anomalies simultaneously for any of the models considered. Models VII Sj0 and Si do, however, lead 
to enhanced Ap B , while not modifying the ti differential spectrum. At the same time they give new 
CP violating contributions to B s and mixing of the right order of magnitude to yield the observed 
like-sign dimuon asymmetry. 

The interaction Lagrangians for the color triplet and sextet scalar fields are given in Q. For vector 
fields the Gf symmetric interactions are given by q^T a q'V^, where T a represents a product of generators 
of color, flavor, and weak SU(2), or some subset thereof, while are the ur, <Ir or Ql family triplets, 
as appropriate 1 . To write down the Gf breaking interactions it is useful to (initially) adopt some of the 
formalism of MFV and promote Y\j,d to spurions that formally transform as bi-fundamentals of Gf 

Y V -> V v Yu V%, Y D -> V D Y D V ] Q . (2) 

Here Vu,d,q are elements of SU(3)tj,d,q> respectively. Assuming full MFV breaking of Gf, all interactions 
are then formally invariant under Gf even for nonzero Yukawa couplings. We will mostly work to the first 
nontrivial order in top Yukawa insertions (the resummation to all orders can be done using a nonlinear 
representation of Gf, see |[T3l[T9] - l2ll ). The explicit forms of the interactions for all the vector models are 
given in Appendix |A| Here we show two examples, models VT Sj0 and Xg 6 . 

Fields VT S ( ) are SU(3)tj r flavor octets, color singlets (octets). The individual field components are V B 
(V A ' B ), where the color label A and flavor label B both run over 1..8. To compress the expressions we 
introduce 

V£ = T B Vf, V° = T A T B V A > B . (3) 

with flavor (color) Gell-Mann matrices T B (T A ) normalized to Tr[T A T B ] = 5 AB /2. The renormalizable 
interactions between quarks and fields VI S)0 are 

Cyr= 4'° u R f s '°u R . (4) 



2 In Hp symmetric models the Lagrangians are trivially obtained from the corresponding Gp symmetric Lagrangians, allowing 
for the possibility that only particular submultiplets of the Gf symmetric representations are present. 
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There are also terms that break Gf — > Hp, 

A£vi s , = [r?2° u R (f s '°^u)u R + h.c] + rjl'° u R {A v f s ' Au)u R + .... (5) 

We kept only the breaking due to Ajj = Y v Yjj insertions. These are diagonal in the up-quark basis and 
change the coupling of third generation quarks to the vector fields (explicitly written out in Appendix |A| 
Eqs. ( |A5[ ), ( |A6[ )). Note that 772'° can be complex and possibly a source of beyond the SM CP violation, of 
interest when considering B s mixing, while fj^ is real. Insertions of Y^Y^Y^Yy are also possible to break 
the symmetry further and are almost diagonal in the up quark basis, while the off-diagonal elements lead to 



FCNC's. We postpone the discussion of these until Section IV D 

The Gf — > Hp breaking also splits the vector mass spectrum. The flavor invariant mass terms are 

C^fJ = (1 + 6 a>0 ) \m 2 Va Tr V£ V^ a + A (H^ H) Tr V^ a }, a = o,s, (6) 

where the color and SU(3) indices are suppressed, and there is no summation over a = o,s (the Kro- 
necker delta 5 a)Q insures the proper normalization for the color octet fields). Note that V is defined 
when rotating to the mass eigenstate basis; see Appendix A. Adding the Gf — > Hp breaking terms 
(1 + 5 a ,o)(imlrTr[V ll A u V' 1 ] and (1 + 5 0)0 )C 2 m^Tr[A[/ A V V^] , the mass spectrum of the vector 
states is (suppressing the o, s labels) 

o 2^2 2 22C12 2 22 4 m\ 

™i,2,3 = m v + -v , m 4)5i6i7 = m 1 + m v — y t , m 8 = m 4 + m v — y t . (7) 

Note that A, (1,2 are all real. The vectors Vj.,2,3 and 14,5,6,7 are degenerate since SU(2)u is only broken by 
light quark Yukawas, not by yt- 

Fields Xg 6 are weak doublets in the bi-fundamental representation (1, 3, 3) of the flavor group. The 
color anti-triplets have field components (V^)^, and color sextets the field components (1^t)[ JQ( g = 
(Vfj,)ij p a , with r the weak SU(2)l index, a, /3, 7 the color indices, while i and j are the indices of the 
(1, 3, 1) and (1, 1, 3) representations respectively. The tree level quark coupling Lagrangian terms are (sup- 
pressing all the indices apart from color, see also Eqs. ( A30 ), ( |A31| >) 



C Xs = vi e af)y 4F 7 Ql + h.c, C Xe = m Qf + h.c. (8) 

Note that the fields transform as — > VdV^Vq, where cIr — > VdcIr and Q C L — > VqQ c l . The mass 
terms are 

C™£ = mlTr(V„V^) + A (fft H^V^V^) + X'Tv{H r V^H^V^) + X'Tr(H^V^H Cs ^V^), (9) 

where we have suppressed all the traced over flavor, color and weak indices (except in the last two terms 
where we show explicitly the weak contractions). We use [(V^)^]* = (V/l)r7, and similarly for the sextet. 
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Note that the last two terms break the mass degeneracy between charge 1/3 and charge —2/3 components 
of Vfj, weak doublets. The flavor is broken through Yukawa insertions 

AC^/m 2 v = CilV(V^Y D \M) + C 2 Tt{Y d Y ] d V^) + ^{V^Y^) + • • • , (10) 

where we do not write down the terms with more than two Yukawa insertions or the similar terms with 
Higgs fields. The resulting Hp symmetric spectrum for charge 2/3 vectors is 

= "* + + A > 2 - = ™" + <U) 

The interactions of mass eigenstates Vki with mass eigenstate quarks (denoted with primes) are given by 
(showing explicitly only color contractions, see also Eqs. ( |A32| ), ( |A32[ )) 



Cx- S = m (4)°fi 7 ( u 'l) cP + m £*m {d' R ) a ?2VcKM {& + h.c, 

£x 6 = vi (4) a (f W K) c/3 + m (d 7 R r(f) a ^v C Ku K) c/3 + h.c, 

where V% and V2 are the mass eigenstate vector fields of the SU(2) doublet. The residual Hp flavor uni- 
versality of these interactions can be broken by insertions of the spurions Yy Y v and Yp Y D . In MFV this 
is the only form of further flavor breaking. The rest of the Lagrangian constructions are collected in the 
Appendix [A] 



III. PHENOMENOLOGY OF TEVATRON ANOMALIES 



We now discuss two recent experimental anomalies observed at the Tevatron: the large forward- 
backward asymmetry Ap B , and the like-sign dimuon anomaly in B s decays. In this section, we systemati- 
cally address the following questions: 

• Is it possible to explain either of the two anomalies assuming Hp symmetric models? By which 
charge and flavor assignments? 

• Are closely related experimental constraints simultaneously obeyed? 

• Is it possible to explain both anomalies using just a single Hp symmetric field? 

A common feature of models put forward to explain the Ap B anomaly [£3l [T5I4T81 122| - |831 is that they 
have 0(1) NP couplings to the up quark. They fall roughly into two classes, those with s-channel exchange 
above a TeV 3 ||25j|26l|75l, in which case the axial vector NP couplings to the top quarks and up quarks must 



3 For a recent exception with a sub-TeV axigluon, see 1171 . 
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be of opposite sign |H3|26l, and those with sub-TeV t-channel exchange |[T8ll22ll23ll26ti28ll57l . in which 
case large inter-generational couplings are required (for additional possibilities, see [29]). The couplings in 
either class could arise from large flavor violation in the underlying theory, which may lead to violations of 
FCNC constraints in K° — K°, B^ — B®, D° — D° mixing, B — > Kir, or b — > 57, unless the couplings 
are carefully aligned (see, e.g., l30l 13111341 ). Moreover, the i-channel models can lead to excessive (flavor 
violating) single top or same sign top pair production at the Tevatron and LHC [22]. 

However, flavor violation is not necessary for large A l p B 11351 . Note that in tt production no net top quark 
flavor charge is generated. Furthermore, models with an unbroken Hp subgroup do not generate FCNCs 
in processes with light quarks. The exact size of FCNCs then depends on the size of Hp breaking. If this 
breaking is MFV-like the FCNCs are generically suppressed below present bounds. The flavor symmetries 
also eliminate single top and same sign top production. 

Many new models have also been put forward to explain the D^dimuon anomaly ll84l - [9~T1l . Together 
with possible indications for deviations from the SM in B s — > J /4> 4> decays and B — > t~ v decays, it 
may point to a NP phase in B s ^ mixing. Intriguingly, MFV suffices to explain the dimuon anomaly [92]. 
After discussing flavor symmetric fields and the A l p B anomaly, we will examine whether these fields can 
also give large enough contributions to B s mixing, under the assumption of MFV breaking of Hp. 

A. General analysis of the t i forward backward asymmetry 

We entertain the possibility that A l p B is enhanced above SM levels via tree level exchanges of fla- 
vor symmetric scalars or vectors. The experimental evidence for such enhancement is as follows. Using 
5.3fb _1 of data CDF measured an inclusive asymmetry A§ B = 0.158 ± 0.072 ± 0.017 in the tt rest 
frame (fixing m t = 172.5 GeV) (SJ. In a channel with both t and t decaying semileptonically an even 
larger asymmetry was found, A FB = 0.42 db 0.15 db 0.05 H. Similarly, a recent D0 analysis finds 
A FB = 0.196 dz 0.060ltoo26 usm g 5.4fb _1 of data Q. Combining in quadrature the statistical and sys- 
tematic errors of the three measurements gives A l p B = 0.200 dz 0.047. This is to be compared to the SM 
prediction A F ' B M = 0.072to:oo7 from an approximate NNLO QCD calculation [§3 with m t = 173.1 GeV 
and using the MSTW2008 set of PDFs [94]. Inclusion of electroweak corrections leads to an enhancement 
of the asymmetry, with A l p B M = 0.09 ± 0.01 recently obtained in Ref. [95]. In the pp frame, a recent 
approximate NNLO calculation (96J gives A§ B M = 0.052+ ) ; ) ) g with m t = 173 GeV, to be compared 
with the CDF value of A§ B = 0.150 ± 0.058 d= 0.024 10. The approximate NNLO SM predictions use 
the known NLO results ||97ti99l and build on recent progress in NNLO calculations II100H1051I . D0 also 
reports a leptonic asymmetry A l FB = 0.152 dz 0.038^q'q{3 to be compared to the MC@NLO prediction of 
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Observable 


Measurement 


SM predict. 




0.158 ± 0.072 ±0.017 E| 






Att 
FB 


0.42 ±0.15 ±O.O5|0 


> ~ 0.200 ± 0.047 (7.24t£$±g;;$) • 10~ 2 ||93l 




0.196 ± 0.0601^26 Q 






Af B (M a > 450GeV) 


0.475 ± 0.101 ± 0.049 




A§ B (M t t < 450GeV) 


-0.116 ± 0.146 ± 0.047 


(s-2±8:S) • 1(r2 123 


A%(\Ay\ < 1.0) 


0.026 ±0.104 ± 0.056 @ 


(4-77i°j9) • lO" 2 El 


A%(\Ay\ > 1.0) 


0.611 ±0.210 ±0.147 


(14.59t 2 ;^) • lO" 2 El 




(6.9 ± 1.0)pb fl08l 


f (6.63i°;°° 7 )pb 1106) 
\(7.08^;°m 6 7 ) P bll07j 



TABLE III: Measurements and predictions for observables in tt production at the Tevatron. We quote the approximate 
NNLO QCD prediction of A^ B from ||9l using MSTW2008 PDFs J94). The two other choices for PDFs give results 
in agreement with these [93]. Among the cross section predictions obtained in 1 106 1 we quote the IPIscet one. 



A 



,SM 
FB 



0.021 ±0.001 Q. 



CDF reported evidence that the anomalously large asymmetry rises with the invariant mass of the t 1 
system, with A$ B {M$ > 450 GeV) = 0.475 ± 0.114, while Af B (M$ < 450 GeV) = -0.116 ± 
0.153 [6]. A similar rise of the asymmetry with respect to the absolute top vs. anti-top rapidity difference 
\Ay\ = \y t - yi\ was also reported by CDF with Af B (\Ay\ < 1.0) = 0.026 ± 0.104 ± 0.056 and 
A§ B (\Ay\ > 1.0) = 0.611 ± 0.210 ± 0.147 0. The recent D0 analysis Q does not observe a significant 
rise of the "folded" detector level asymmetry with respect to M 4 j and | Ay\. However, until these results are 
unfolded they can not be directly compared to the CDF measurements, although at the detector level they 



appear to be consistent within errors. We collect the above results in Table III 

Any NP enhancement of A l p B must not spoil the agreement between the measured production cross 
section, a t j, and the SM predictions. At NLO with NNLL summation of threshold logarithms, the SM 
prediction is a a = (6.63±Q™)pb lfl06ll (using MSTW2008 pdf sets and IPIscet choice of kinematic 
variables and resummations - the other choices give consistent results but with larger error bars). This is 
somewhat smaller than the approximate NNLO result (for m t = 173 GeV), a t i = f.OS+o^to^r P b 1110711 
(see also [lOOl 11041 11051 ). Both of these results agree well, within errors, with the measured CDF result 
based on 4.6/6- 1 of data |[T08l a t - t (m t = 173.1 GeV) = 6.9 ± l.Opb. Thus the NP contribution to the tt 
cross section, crf^, is tightly constrained. 

Good agreement between experiment and SM predictions is also seen in the differential cross section 
da/dM t f. This has important implications for the viability of different NP models. For instance, comparing 
the measured and predicted cross sections together with the measured and predicted A f BB , for M t i > 450 
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GeV, one finds that the NP contributions need to reduce the backward it-scattering cross section (a statement 
valid at 2a). This can only happen if NP interferes with the SM ll35ll . NP in the s-channel which interferes 
with the single gluon exchange amplitude must therefore be due to color octet fields. In general s-channel 
resonances lead to significant effects in da/dM t ^. However, this may be avoided for a purely axial gluon 
that is broad fl7l . in particular regions of parameter space. There are no such clear requirements on the 
charge assignments of possible t-channel NP contributions. However, a characteristic high mass tail in the 
spectrum could lead to tension with the Tevatron and future LHC cross section measurements at large M t j. 

We collect expressions to be used in our analysis below. The total cross section a t j, forward-backward 
asymmetry A FB , and the cross sections for forward and backward scattering op,B are defined as 

At f of — &b f 1 ^ , n da 

Oti = OF + OB, A FB = , a F{B) = dcos9- (13) 

a F + (TB Vo(-i) acosfl 

where is the angle between incoming proton and outgoing top quark. We use NLO SM predictions for a t j 
and do t ij dM t i, and LO predictions for the NP corrections (including interference with the SM). To obtain 
A l FB we define a partonic level asymmetry, 

aNP+SM g> -B , ASM ( ° \ (U \ 

^FB /^.NP+SMn , / NP+SM\ fB \ rrSM i n NP ) ' V ' 

{(T F )LO + (^B >LO \ a + a J 

which is to be compared against the binned unfolded partonic level results of [6]. We use the NLO + NNLL 
SM predictions for the forward, backward and total cross sections 11091 , the ti spectrum H109II and A f p B M 
1931 . For concreteness, for M t - t < 450 GeV we take the central values a SM = 4.23 pb and Ap' B M = 0.052, 
while for M tl > 450 GeV we take a SM = 2.40 pb and A%' B M = 0.111 (MSTW08 pdfs); for the inclusive 
asymmetry (in the ti rest frame) we take Ap B = 0.0724. 



1. Acceptance effects 

As pointed out in lT30l 1 1 1 01 . care is needed when comparing NP predictions to the experimental parton 
level A FB and M t i differential spectrum deduced by CDF [6, 108], since the deconvolution was done 
assuming the SM. The acceptance corrections are especially important if NP enhances top production in 
the very forward region This is because the SM ti event distribution is more central. We take into account 
the CDF experimental cuts using correction factors e^. For the i— th bin in M t i one needs to multiply the 
calculated partonic ti cross section da^ p / dM t [ by 6{ in order to compare with the CDF measured partonic 
cross section (da NP / dM t t) CDF 

'da NP \ CDF fda" p \ 



dM H 



£ - x Kan ■ <15) 
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M«(GeV): /Ay: 


0-0.6 0.6-1.2 


1.2- 1.8 1. 


8-2.4 2.4-3.0 


350 - 400 


2.42 


2.23 


- 


_ 


_ 


400 - 450 


3.44 


2.74 


1.95 


— 


— 


450 - 500 


4.29 


3.32 


1.75 


- 


— 


500 - 550 


5.25 


3.50 


2.34 


1.18 




550 - 600 


5.61 


4.39 


2.67 


1.15 




600 - 700 


6.59 


4.81 


3.01 


1.09 


0.42 


700 - 800 


7.38 


6.06 


3.49 


1.34 


0.45 


800 - 1400 


6.68 


6.26 


4.30 


1.83 


0.42 



TABLE IV: Acceptances (in %) to be used for emulating the CDF deconvolution to parton level measurements (cf. 
eqs. ( fl~5] l - (fTT)). Note that the do not depend on the NP model and are symmetric in Ay O —Ay. 

There is no summation over i in this equation. Since CDF is using SM acceptances and no bins in Ay in 
the deconvolution of the da/dM t i measurement in 0, the are given by the ratio of acceptances for the 
NP model and the SM 

,-NP 



(16) 



where e ^ p ( SM ) are calculated by splitting each i-th M t j bin into j bins in Ay = yt — Hi 

V t NP(SM) 



NP(SM) 



NP(SM) ' 



(17) 



and the sum is over the bins in Ay. Here ejj is the acceptance for each {M t i, Ay) bin, and (j^ p ( SM ) [ s th e 
corresponding cross section integrated over the bin. The above expressions are approximate in so far as 
the bins have finite sizes, and the spill-over of events between different bins is not taken into account. The 
acceptances are calculated by simulating the partonic t — i cross section using MadGraph4 .4.30 111 1 IB . 
decaying the top quarks in Pythia6 . 4 H112H . which also simulates the LO showering and hadronization, 
and using PGS for detector simulation. The events were read into Mathematica, where the cuts from 
(6) 1 1 1 311 were implemented. The resulting values for the acceptances are collected in Table IV 



In Fig.[T]the correction factors ej are shown as a function of the M t j bin for two benchmarks points (cor- 
responding to illustrative couplings and mass valuers) in models VI S)0 which exhibit substantial departures 
from the SM acceptances. Results for da / dM t j and A^ B for the two model VI Q benchmark points in Fig. [I] 
were shown previously without acceptance corrections 051 . The differential distributions with acceptance 
corrections are compared to those in ||35ll in Fig. [2] We see that the corrections bring the predicted spectrum 
for the light vector example into good agreement with experiment. For the remaining Gf representations 
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FIG. 1: The correction factors q for i-channel dominated vector models VI S (left) and VI D (right), where 
(m v , v^fe?733, Via, r y /m y )=(300 GeV, 1, 1.33, 0.08)[dashed red]; (1200 GeV, 2.2, 4.88, 0.5)[solid blue]. 

discussed below, the correction factors are not as important. For models Sy, VIII S and IX S the e» are below 
15%, and for the others they are below 5%, for all M t j bins and all benchmarks points considered. 

The pattern of acceptance corrections can be understood from the angular dependence of the NP 
contribution to the differential cross section in each model. For instance, the t-channel exchange of a 
vector with mass my leads to a Rutherford scattering peak in the forward direction for m?r >> rriy. 
Specifically, the expressions for the NP cross sections contain characteristic t-channel (1 — cos 6>) 2 ' 4 factors 
in the denominators, whose angular dependence is reinforced by (1 + cos6>) 2 factors in the numerators, 
where 6 is the top quark scattering angle in the tt center of mass frame. Thus, models VI S)0 with light 
vector masses favor forward top-quark production at large M t j, yielding t\ that are substantially less than 1 
in the high M t i bins. In models Si (t-channel) and Sy, Syi (n-channel), the angular dependence introduced 
by the characteristic (1 =F cos#) 2,4 factors in the denominators is offset by (1 =F cos6>) 2 factors in the 
numerators, which leads to central NP top-quark production, as in the SM. The result is that the in these 
models are actually slightly larger than 1 . 

To apply the acceptance corrections to A FB we follow CDF, where A l p B was obtained using four bins 
in M t i and Ay: M t j above or below 450 GeV and Ay positive or negative. The correction factor for each 



of the four bins is given as in Eqs. (p3]>-(^7|>, except that the sum in ( [17] ) now runs over all j with either 
Ay > or Ay < 0, and over the appropriate values of i with either M t i > 450 GeV or M t i < 450 GeV. 
We find that the corrections are small for all of the models we consider. For instance, for the light vector 
color octet example in Fig. [I] the shift is from an uncorrected A l ^ B = (0.10949,0.357) to a corrected 
A l p B = (0.10953, 0.339), where the first and second numbers are the low and high mass bins in Ap B . The 
small shifts in Ap B are due to the coarse binning in M t i II 101 . In particular, the high mass bin is dominated 
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FIG. 2: The effect of acceptance corrections on the vector model discussed in 051 . Before acceptance cor- 
rections (left), after corrections applied (right). For the parameters (my, y/VijV33i Va^v l m v)'- solid red 
(300 GeV, 1, 1.33, 0.08); dashed blue (1200 GeV, 2.2, 4.88, 0.5). 



by events with il% near 450 GeV, which are more central, as in the SM. 



2. The 1 1 phenomenology of Hp symmetric scalar fields 

The flavor symmetric models introduced in Section [n] and collected in Tables [TJ |n] can couple to light 
and heavy quark bilinears with unsuppressed couplings. They are thus interesting candidates to explain the 
A l p B anomaly, as noted in |[3j[35]|36l. In the case of scalars, SU(2)l singlet color triplets or color sextets 
and SU(2)l doublet color singlets have previously been identified as being promising for explaining the 
A FB anomaly |[T8l[32l[33ll35]437llll4ll . Here, we will focus on some of the flavored versions listed in Table 
PH Our results overlap with past studies, but we also include color triplet and sextet scalars that couple to 
initial state down quarks, that have not been studied as extensively. We find that these models may also be 
viable, although they generically require a coupling that is a factor of ~ 2 larger than if the up quark is in 
the initial state. This rule of thumb also holds for the vector models that we will study in next subsection. 

We show results below for models Sv,vi that couple to initial state up quarks, for models Six,x that 
couple to initial state down quarks, and model Si that couples to initial state up and down quarks. (Flavored 
color sextets and color triplets were considered previously in |[35l |36l ). We expect models Sxi — Sxiv 
to yield results similar to those for Six,Sx- The scalar models Sm and Srv are in the 1 and 8 color 
representations respectively. These models are known to generically suppress A l p B when interfering with 
the SM |[T8l and we do not consider them further. 



17 



-a 
a 



10 

9 
8 
7 
6 
5 



Scalar field V 



Scalar field V 



\ 



200 400 600 800 1000 1200 1400 

M s [GeV] 

Scalar field VI 



-a 



10 

9 
8 
7 
6 
5 
4 



1 



-5 



100 
10 

1 

0.1 
0.01 



200 400 600 800 1000 1200 1400 

M s [GeV] 

Scalar V,VI 



Scalar V. 
Scalar VI - 



400 600 800 1000 1200 1400 
M tt [GeV] 



0.4 



0.3 



X 0.2 



0.1 
0.0 

0.4 
0.3 



^ 0.2 
0.1 



0.0 



/ \ 



* 



200 400 600 800 1000 1200 1400 

M s [GeV] 

Scalar field VI 



/ 
i 







IL 

200 400 600 800 1000 1200 1400 
M s [GeV] 

Scalar V, VI 



0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 



Scalar V 
Scalar" VI " 



200 400 600 
M tt [GeV] 



800 



FIG. 3: Upper two rows: predictions for the inclusive cross sections a(ti) and inclusive forward-backward asymmetry 
ApB(tt) as a function of scalar mass mj for models Sy, Syi and couplings rj — 1/4 (solid line), l/2\/2 (dotted), 
1/V2 (dot-dashed), 1 (dashed) compared to la and 2a experimental (shaded) bands. Predictions for A l p B with M t i 
below and above 450 GeV and for da t i/dM t j are compared to the experimental data from Il6l ll08 l in the last row for 
benchmark points labeled with a * in the inclusive predictions (ms = 400 GeV for Sy and 800 GeV for Syi). 



The interaction Lagrangians for scalar models Si and Six are 

Ci = r] ( S^Ui L uj R + S'jdiLUjR) + h.c, 
Ax = rie^Sfjd 1 ^ + /i.e., 



(18) 
(19) 
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where i,j and a, /3, 7 denote flavor and color indices, respectively. Explicit forms of the interaction La- 
grangians for the remaining scalar models under discussion, Sv,vi,X w& collected in Appendix [A] In the 
Gf symmetric limit the decay widths of the scalars are nrf ms /167T, where r\ is the coupling of the scalar 
field to the SM quarks, and k = 1, 8, 8, 2, 2 in models Si, Sy, Syi, Six, Sx, respectively (assuming all 
quark decay channels are open and ignoring phase space effects). The interaction Lagrangians explicitly 
defining rj are collected in Appendix [Aj while the relevant NP cross section formulae for top quark pair 
production can be found in Appendix B. 

In Figs. [3] and [4] we collect predictions for the inclusive A^ B and a t 1 as functions of the scalar masses 
in models Sy, Syi, Six, an d Sx, for several values of the couplings 77 to quarks. These results do not 
depend on whether one is considering the Gf symmetric limits in Eqs. ( fT9] >, ( |A37 1, and ( A38 1), or the 
Hp symmetric limits (with the couplings of the light quarks to the top identified with the chosen values of 
rj). Large effects on A BB are possible, while remaining within the 2a bounds for a t t- In general 0(1) 
couplings rj are. required. Therefore, some Gf — > Hp flavor breaking will be necessary in order to satisfy 
the dijet constraints on the couplings to light quark pairs, particularly in models Syi and Sx (we comment 
on this further below). For example, recent LHC dijet measurements imply that in model Syi the couplings 



of light quarks to ~ 1 TeV mass scalars should be < 0.1, see Section IV C 



A strong constraint on these models is consistency with the measured da{tt) / dM t i distribution. To 
demonstrate the potential of these models to explain the Ap B anomaly and their impact on this spectrum 
we pick particular benchmark values for the scalar masses and couplings, denoted by a ★ in Figs. [3] [4] 
The resulting da(tt)/dM t j spectra and A FB in the low and high il% bins are shown in the last rows of 
Figs. [3] [4] We do not apply K-factors to the NP cross sections in any of the models we consider. A 
genuine NLO calculation for these models is beyond the scope of this work, but is clearly required before 
precise conclusions can be drawn regarding the effects of the various models on the do jdM t i distributions. 
All chosen points predict a significantly enhanced A l p B in the high invariant mass bin. In general, these 
models reduce the good agreement between the SM and measured da{tt)j dM t i spectrum, while improving 
the agreement in A FB . Clearly, there is significant tension between satisfying the constraints on the M t i 
distribution while simultaneously creating a large A l FB enhancement, as has recently been emphasized in 
|[37l l44l . Recent D0 data [7] shows a preference for a smaller A FB in the high M t i bin. While the M t i 
dependence of A FB has not been unfolded by D0 , and thus cannot be directly compared to the results in 
Figs. |3]and|4| some of the tension may be alleviated with the new data. 

The flavor breaking naturally present in this framework can reduce the remaining tension with measure- 
ments associated with the light quark sector, e.g., from the dijet bounds mentioned above. For instance, 
keeping the leading insertions in that break Gf — > Hp, the Lagrangian for scalar model Six in the mass 
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FIG. 4: The same as Fig. [3] but for models Sx,ix and couplings -q = l/\/2 (solid line), 1 (dotted), \/2 (dot-dashed), 
2 (dashed). Benchmark masses are mj = 750 GeV (Six ) and mg = 800 GeV (Sx )• 

eigenstate basis becomes, 



c a = m {d' R ) ai {u' R ) pj + { m + 2 m yl) (d' R ) ai (t' R )p s'^ 



(20) 



with j = {1, 2} and i = {1, 2, 3}. The mass spectrum is split, with mf 3 = mf- + Ci77t|. Depending on the 
signs of 771,2 and (1 one can have an enhanced coupling to top quarks (compared to couplings to light quarks) 
with a suppressed mass scale of the mf 3 components of the scalar. This can naturally lead to larger effects 
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FIG. 5: Upper row: predictions for the cross section a(tt) and inclusive forward-backward asymmetry ApB(tt) as 
a function of scalar mass ms for model Si and coupling 77 = 1/3 (solid curve), 2/3 (dotted), 3/4 (dot-dashed), 1 
(dashed). Predictions for A l p B with M t f below and above 450 GeV and for da t t/dM t 1 are shown in the last row for 
the benchmark point labeled with a * in the inclusive predictions (7715 = 110 GeV). 



on top phenomenology, such as A l p B and cr t j, while at the same time weakening the impact of the dijet 



constraints (see Section IV C ). This situation is similarly realized in the other scalar models. Alternatively, 
one could consider minimal Hp symmetric versions of models Six and Sx, which only contain the scalar 
Gf submultiplets transforming as (1, 2, 1) under Hp (they only couple the light quarks directly to the top 
quark), while in the minimal Hp symmetric versions of models Sy and Syi one would consider the (2, 1, 1) 
representations of Hp. 

Finally, we discuss model Si. 4 (Models with flavored scalar doublets have previously been considered 
in |[32l l33l ). For illustration, we can define the left-handed (LH) quark fields in Eq. ( p~8] ) in either the 
up or down quark mass eigenstate basis (with the right-handed (RH) up quarks in their mass eigenstate 
basis). In both cases, the recent D0 upper bounds HI 1611 on anomalous resonant dijet production in W + 



4 Note that two studies have recently concluded 1371 1 1141 that a (color singlet) SU(2)l doublet would be the most promising, 
among the possible scalar representations, for explaining the ti asymmetry with minimal distortion of the Tevatron and LHC 
M t i spectra. 
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jets imply that the couplings of the scalars to light quarks should satisfy rjij < 0.2, i, j = {1, 2}, given the 



light scalar doublet masses favored by the Ap B anomaly (see below and Section IV C I. Significant breaking 
of Gf — > Hp would, therefore, be necessary in order to accommodate O(l) scalar-top quark-light quark 
couplings. However, we reiterate that Hp symmetry protects against NP contributions to K — K and D — D 
mixing, as well as to single top and same sign top production. 



If the LH quarks in Eq. 18 are defined in the up mass basis, then the B — > Kir branching ratio constraints 



071 II 1511 would also impose the bound 773^ <0.2, i = {1,2}, see Section IVD Therefore, in this case 
top phenomenology would be well approximated by a minimal Hp symmetric model consisting of scalar 
SU(2)l doublets which transforms as (1, 1, 2, a) under Hp, where the U (1)3 charge a is opposite to that of 
the top and bottom quarks (U(2)q l is defined with respect to the LH up quark mass eigenstate basis). We 
can write the corresponding interaction Lagrangian, in the quark mass eigenstate bases, as 

= V (5? 3 %l t' R + % (^ckm)M'l4) + h.c, (21) 

where i = {1, 2}, j = {1, 2, 3} and Vckm is the Cabbibo-Kobayashi-Maskawa (CKM) matrix. If the LH 
quarks in Eq. ( p~8] ) are defined in the down quark mass eigenstate basis, there are two minimal Hp symmetric 
and flavor safe alternatives for Ap B enhancement: exchange of scalar SU(2)l doublets which transform as 
(2, 1, 1, a), or as (1, 1, 2, a) ( U(2)q l is now defined with respect to the LH down quark mass eigenstate 
basis). The relevant interaction Lagrangians would be 

= V (S^ u ]L {V CKM ) ri u' lR + b' L u' iR ) + h.c, (22) 



or 



= V (5? 3 rf,L (YcKM^it'n + S7 3 d' iL 4) + h.c, (23) 

respectively, where i = {1, 2}, j = {1, 2, 3}. For the light scalar masses ( < 130 GeV) and 0(1) couplings 
7] favored by A FB , options 1(a) and 1(c) above (Eqs. pT) , ( |23] >) are disfavored at the 4<r level by consttaints 
on non-oblique corrections to the Z couplings to quarks, whereas option 1(b) (Eq. p2])) is consistent at 



2a, see Section pVB In a generalization of option 1(b) to complete Gf multiplets, non-oblique correctons 
would dictate that the couplings 77^3 need to be roughly a factor of 2 smaller than 7731 for i = {1,2}, while 
Wjj constraints would dictate the bounds rjij < 0.2 for i, j = {1, 2}. 

In Fig. [5] we show predictions for A^ B (inclusive and Mg bins) and da / dM t i corresponding to Eq. (22). 
Our results confirm the findings in ll37ll for Ap B enhancement with scalar doublets: given CDF data the 
preferred scalar mass is small, below 130 GeV, while the coupling 77 is O(l). Also, the difference between 
the asymmetries in the low and high mass bins is not large when the latter is enhanced (which could be 
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FIG. 6: Upper two rows: predictions for total cross sections <j(tt) and inclusive forward-backward asymmetry 
ApB(tt) as a function of vector mass my for models II so and couplings f q f t = 1/256 (solid curve), 1/128 (dotted), 
1/32 (dot-dashed), 1/8 (dashed). The bands are the one and two a measurements of Ap B and cr t j. The predictions 
for Apb with M t i below and above 450 GeV and for da t i/ dM t t are shown in the last row for benchmark points 



labeled with a * in the inclusive predictions (My = 500 GeV). The data are summarized in Table III 



welcome in view of the recent D0 A l p B measurements), and there is minimal impact on the M t j spectrum. 
Finally, we note that for our benchmark points the bounds on the top quark decay width are not violated by 
new top quark decays to a scalar and a light quark (even though there are 2 times as many such modes in 
1(b) than previously considered in ll37l ). 
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The 1 1 phenomenology of Hp symmetric vector fields 



In flavor symmetric scalar models, ti production can only proceed via the u or t channels. Flavor 
symmetric vector models differ in this respect, as ti production can also proceed in the s channel or both in 
s and t channels simultaneously. Top phenomenology in these models is dictated to a large extent by which 
channel dominates. The sizes of the couplings required in order to explain the top asymmetry are roughly 
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fixed by whether the vectors couple to up or down quarks in the initial states, as in the scalar models. The 
Lagrangians for the vector models are given in Section [n] and Appendix [a] Expressions for the qq — >■ ti 
cross sections for models I — XI are collected in Appendix [B| They are given in terms of effective couplings 
of the vectors to light quarks (f q ), to top quarks (ft), and to top and light quarks (f q t), defined in terms of 



771^2 etc, see Table VII In this section we use these couplings as numerical inputs. In the Gf symmetric 
limit (yt = 0), one has f q = ft = fqt- Among the vector models, I s and V SiQ do not contribute to top 
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FIG. 9: Predictions for da/dM t j and Afb for a set of vector models, IV Si0 , VI Sj0 , VIII Sj0 , IX s o fixing my =500 
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production, and I s , II S and III S do not interfere with the SM amplitude at tree-level. 

We begin with an examination of the representative models II S , 0) which only contribute in the s-channel. 
The vector decay widths are given by my (2/ 2 + / 2 )/87r (II S ) and my (2/ 2 + _f 2 )/487r (II ), neglecting 
phase space suppression due to top quarks in the final state. In the first two rows of Fig. [6] we collect 
predictions for the inclusive A l p B and er t j as functions of the vector masses, for several values of the product 
fqft- We have taken vector decay widths appropriate for the Gf symmetric limit: my3f q ft/8ir (II S ) 
and myfqft/16^ (II ). In the last row, predictions for dajdM t i and A FB in the low and high M t i bins 
are shown for the benchmark points (denoted by a Good fits to the inclusive cross section and the 
top asymmetries can be obtained. However, agreement with dajdM t i is much harder to achieve, as the 
effect of the s-channel resonance is clear. The bounds from the Tevatron can be avoided if my > 1 TeV. 
Unfortunately, for my > 1 TeV and O(l) couplings the effect on A FB is strongly reduced. The LHC 
da/dM t i spectrum would also be problematic. It is possible, though, to hide a light s channel resonance 
via an increase of its decay width due to additional non-qq final states, as has recently been discussed in 
the context of pure s-channel axial- vector models lfl7l . Another possibility for increasing the vector decay 
widths may be to consider the Gf breaking hierarchy f q « ft, keeping the product f q ft fixed. 

The flavor octet models VI S and VI Q contain admixtures of s and t channel contributions, while the 
flavor (3, 3, 1) models VII S and VII Q are purely t channel. The s channel contributions in VI SiQ are due 
to the vectors associated with the generator T 8 of SU(3)u R and the coupling product f q ft. The t channel 
contributions are due to the vectors associated with the generators T 4 ' -' 7 and the coupling product fl\. While 
the generic problem of s channel resonances persists, it is mitigated by their relative suppression compared 
to the t channel contributions (the coefficients are significantly smaller than C3 for these models, see 
Table VII| ), by possible cancelations between the two channels, and by the absence of interference between 
the SM and NP s channel contribution in VI S (C5 = 0). The s channel effects can be further diminished, 
if the vector decay widths are enhanced by additional non-qq final states. 5 This assumption was implicit in 
the model VI Q examples presented in |[35l and Fig. [T] In the following, we only consider vector decays to 
quark bilinears. In the Gf symmetric limit (f q = ft = f q t) the decay widths can be written as f qt my /167T 
(VI S and VII S ) and /^my /96tt (VI d and VII Q ). Finally, s channel effects in VI Sj0 would be suppressed 
in the Gf breaking limit f q f t « f qt (this breaking is natural due to the large top yukawa), and would 
be absent entirely in the minimal Hp symmetric realizations which only contain the complex SU(2)u R 
doublets associated with the SU ($)u R flavor generators T 4v ' 7 . 



5 Given that these models are by themselves non-renormalizable, it is reasonable to expect that renormalizable UV completions 
could contain additional vector decay modes. 
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We show the effect of the models VI Sj0 and VII Si0 on A l J B and a t i in the Gf symmetric limit in Figs. [7] 
and [8j respectively. In the purely t-channel models VII so , we see that it is possible to enhance A l J B , while 
at the same time not introducing any visible deviations in the da / dm t i differential spectrum. For models 
VI Si0 , on the other hand, we see that it is difficult to enhance A l p B sufficiently, without obtaining excessive 
s channel peaks in the da/dM t j spectra, as in the model II s o examples above. However, as already noted, 
much better agreement with the measured spectrum can be achieved in this case, if Gf is broken down to its 
subgroup Hp 6 - In Fig. |9]we show such examples, with f q f t « f qt , for various vector models. The decay 
widths of the different vectors in a Gf multiplet will not be equal in this limit. For simplicity, we have iden- 
tified all the widths in each multiplet with the quantities m v f g f t /16irx (3, 1/2, 1,1/6, 2, 1/3, 1/2, 1/12) 
in models IV S , IV Q , VI S , VI Q , VIII S , VIII , IX S , IX Q , respectively (which would hold in the Gf symmetric 
limit, neglecting phase space differences). The rationale for this approximation is that the s channel res- 
onance widths are associated with the coupling products f q and / t , whereas the widths of the t channel 
resonances, controlled by the product fl\, have minimal impact on the M t j spectrum. In addition, as we 
will see in Section IV C (see Table VI ), dijet constraints on the light quark vector couplings f q are relatively 
mild, so that the relation fi 1 = ff = f q f t would be phenomenologically viable in the light vector examples 
of Fig. 



4. Limits from LHC measurements of the 1 1 invariant mass spectrum. 

Recently, the first measurements of the tt invariant mass spectrum at the LHC have been presented 
by the ATLAS and CMS collaborations. With up to 0.7 pb _1 the ATLAS collaboration obtains a t i = 
176 ± 5(stat)!jo(syst.) ± 7(lumi.) pb, while using 36pb _1 the CMS collaboration obtains a combined 
cross section o~ t i = 158 ± 10(uncor.) ± 15(cor.) ± 6(lumi.) pb 1112011 . This is to be compared with the 
NLO+NNLL predictions or the approximate NNLO predictions for the inclusive cross section that are listed 
in Ref. 11061 . and range between a t i = 157± 13 pb and a-g = 145± 11 pb, depending on the approximation 
used. ATLAS and CMS have also presented promising first measurements of the differential cross section 
at 200 pb" 1 and 886 pb" 1 [ 121 ], respectively. 

In the near future the differential cross section measurements will place meaningful constraints on the 
models we have considered. We therefore present below some representative examples of their impact on 

6 Whereas in models VII S , the Gf — > Hf breaking is not needed in order to obtain good agreement with the da / dm t f spectrum, 
it may be required to evade dijet and B — > Kn constraints. 
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the differential spectrum. In Fig. 10 we show the ratio 

_ da(tt) N p +S M/dM a 
u ~ da{tt) SM /dM tl 

in models II , VI Q , Sy, and Syi for the benchmark points shown in Figs. [6j [7] and[9j In the ratio we use 
SM NLO+NNLL prediction from Ref. 1110911 , while NP contribution including interference with the SM is 
calculated at LO, mirroring our procedure for Tevatron predictions used in previous sections. In Fig. 10 a 
sharp resonant peak is clearly visible for the disfavored pure ,s-channel model II Q . Once convoluted with 
the experimental resolution the peak will be less prominent. For instance assuming 10 GeV resolution the 
peak leads to an O(30%) enhancement of the cross section in the bin containing the peak. The model VI 
example contains a rise in the tail region associated with the Rutherford scattering peak that is characteristic 
of ^channel models with 0(1) couplings. The scalar models Sy and Syi, being n-channel models, display 
a relative enhancement throughout the entire differential spectrum B4ll . similarly to the case at Tevatron (see 
Fig. [3]). Model Si displays virtually no deviation from the SM spectrum, as pointed out in Refs. ||37ll68l . 



B. Contributions to B s — B s mixing 

There is some evidence for NP contributions to the measurement of the B s mixing phase. A 3.9a 
deviation from the negligible SM prediction has been measured in the like-sign dimuon charge asymmetry 
by the D0 collaboration (HOH. This result is in agreement ll92~lll22H with a hint for a nonzero weak phase 
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in B s mixing (measured through flavor tagged decays II1231I12410 . The hints for NP in B s mixing have two 
preferred solutions ll92l . with 

h s ~ 0.5, a s ~ 130° and h s ~ 2, cr s ~ 100°. (25) 

Here, /i s ^ and 2cr Si d denote the magnitude and phase, respectively, of the NP contribution to the B s ^ mixing 
amplitude, normalized to the SM one (see also ( |C3| >). Note that the above results have been obtained using 
the older measurement of the dimuon asymmetry [9]. Since the two measurements are consistent, we do not 
expect a significant change in the position of the two minima. There is also a slight preference for ha ~ 0.2, 
ad ~ 100°, but at slightly more than la is consistent with zero. At 3a one finds < 0.5 for all a^ [92]. 

The flavor symmetric models can be grouped by whether or not they can give contributions to B s — B s 
mixing. The models in which the new fields only couple to ur fall into the latter category These are models 
II Sj0 , VI Sj0 and Sy.vi- Then there are models in which the NP contributions first arise at loop level: VII Sj0 
and Six,x> which couple to cLr and ur, , XI§ 6 and Syi that couple to ur and Ql- For the remaining models 
tree level contributions to B s mixing are possible. 

For a quantitative analysis one needs to specify the Hp breaking. We will adopt the MFV hypothesis, 
assuming that the breaking in the NP sector is only due to the SM Yukawa couplings, (for an example with 
maximal breaking of Hp see EH). This hypothesis ensures that the FCNCs generated from the exchange 
of NP fields have SM CKM suppression, making the FCNC constraints easier to satisfy. MFV scalar fields 
that can lead to the dimuon anomaly have been discussed in the literature previously [88-91 ]. We begin by 
classifying the flavor symmetric vector and scalar fields that can contribute to B s mixing at tree level. Most 
of the details are relegated to Appendix |C| 

The MFV models fall in two categories:"universal models" leading to class 1 AB = 2 mixing operators 
in the terminology of 021, and "yukawa suppressed" models that give rise to class 2 operators and are 
additionally suppressed by light quark yukawa couplings, y s ^. In universal models the NP contributions 
(normalized to the SM) are equal in the B s and Bd mixing amplitudes. This is in some tension with exper- 
iment, where there are indications for large effects in B s mixing, i.e., 50% — 200% of the SM amplitude. 
Measurements in the Bd system, on the other hand, can accommodate a NP contribution of up to ~ 20%. 
Universal models can thus explain the dimuon anomaly only if the real effect is on the lower end of the 
experimentally preferred range for B s mixing and on the upper range for Bd mixing. 

The simplest examples of universal MFV models that can explain the dimuon anomaly are the two 
scalar doublet model, where in our notation the higgs linear combination without a vev is Sh I89TI911, 
and the color-octet, weak-doublet scalar model Ss flU - Effects are large even for raj = 1 TeV, if the 
couplings of the flavor breaking terms r\{ and the bottom Yukawa coupling yj, are 0(1). For the SM value 



30 



Ub ~ 0.02 the mass scale for the new uncharged scalars would have to be quite low. Such a low mass scale, 
surprisingly, is consistent with present phenomenology for Ss H125I . Similar comments apply to the other 
universal models, III Si0 , IV S)0 , VIII so , K Sj0 and Sxiv- m order to obtain the right size of h s for the dimuon 
anomaly one needs r/j ~ 0.1, if mgy = 1 TeV. 

The yukawa suppressed models can have mixing amplitudes suppressed either by one power or two 
powers of light yukawas. For y 2 s suppressed models, I s . , V Sj0 and Sviii, the B s mixing anomaly can be 
potentially explained only if r\i ~ 0(1) and the masses of NP fields are no more that a few 100 GeV. The 
NP contributions to Bd are predicted to be too small to be observed in this case, since they are additionally 
suppressed by y^/vl- The singly y s suppressed models, Xg 6 and Si, Sni> gi ve contributions to h s of the 
right order of magnitude for r\{ ~ 0(1) and msy ~ 1 TeV. There is also an effect predicted for Bd mixing 
with size hd ~ 0.1. 

A number of flavor symmetric models can potentially explain both the B s mixing and A FB anomalies 
simultaneously. If enhanced B s mixing is due to tree level exchange, the NP fields need to couple to left- 
handed doublets if they are to also contribute to ti production. The relevant models are III S)0 , IV S)0 , VIII s o 
and IX Sj0 , where for TeV masses, the flavor breaking couplings need to be r\i ~ O(0. 1). However, none of 
these models seem to fit the observed data in top pair production particularly well (See Fig. [9] for models 
IV Sj0 , VIII Sj0 and IX s o . The models III s o are pure s-channel and thus lead to large deviations in the ti 
differential cross section.). 

There are also four models of interest that exhibit linear Yukawa suppression and can in principle con- 
tribute to B s mixing and to the A FB anomaly simultaneously. These are models Xg 6 and models Sjii, 
Siv- However, due to their color representation Shi, Siv, have limited potential to explain A l BB , while X§ 6 
reduce rather than increase A l FB . 

This leaves us with models where the contributions to B s mixing are loop suppressed. These are the 
vector models VII Si0 , Xlg 6 and the scalar models Si, Sn, Srx,x- They lead to contributions to B s and Bd 
mixing which are universal, if one assumes MFV Among these, models VII s o and Si can also explain 
the present top data. The contributions to B s mixing can also carry a new weak phase, as required by the 
dimouon charge asymmetry. 

For instance, model Si, see Eq. ( |2~Tj ), contains contributions that are proportional to the product of the 
generally complex quantity 47^773 (where these coefficients are defined in Table |x| in Appendix [c]) with a 
loop function that is of the same order of magnitude as the SM one, with its precise value depending on the 
size of ms- The B — > Kir branching ratios require ^/r/orfe ;$ 1/4 (see Section 
7/0 ~ 1. Thus, CPV contributions to B s mixing of O(20%) are readily attainable for Si scalar masses of 
order 100 GeV. Similar considerations apply to models VII s o . 



IV D 1, while A FB favors 
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IV. EXISTING EXPERIMENTAL CONSTRAINTS ON A FLAVOR SYMMETRIC SECTOR 

We now determine some of the existing experimental constraints on flavor symmetric extensions of the 
SM. We examine bounds on direct production at LEP, from electroweak precision data (EWPD) and from 
collider searches at the Tevatron and LHC. We also determine the (residual) FCNC constraints. 

Our aim is characterize the general phenomenological bounds whenever possible and to check if some 
of the models that we have identified as promising candidates are consistent with other phenomenology 
in some detail. Our results are not completely comprehensive, as at times we have to focus on particular 
sample models. However, we believe our results are valuable to gain some intuition as to the current 
phenomenological constraints on a flavor symmetric sector. 



The LEP constraints depend on how the new states couple to electrons and to the pair of final state 
fermions. In this section we will discuss the constraints on the vector models from LEP in some detail. 
The constraints for scalar models are very similar, see for some details. Limits on anomalous four jet 
events at LEP give a kinematic lower bound for pair production of the scalar models of 105 GeV J3J. For 
the vector and scalar models signals that involve couplings to quarks and leptons are the most problematic 
at LEP. For a Z' vector boson that couples to quarks (electrons) with couplings g z zjB (g z z e A), the bound is 



the order of M z > Z 1 TeV. 

The vectors in nontrivial color or flavor representations (I Q — IV Q or V — XI) are protected from having 
tree level exchanges of this form. For models protected by flavor symmetry one is forced to have the 
insertion of Yukawa matrices. The final operating energy of LEP was \fs = 209 < 2m t so the flavor 
symmetry is sufficient to remove this bound for models V — XI. These bounds do apply to I s — IV S unless 
the couplings to electrons are somehow suppressed, e.g., by extending flavor symmetry to the lepton sector. 

For the models where the bound ( [26] ) does not hold, there are constraints from covariant derivative 
couplings to Z and 7 in the kinetic terms of the vectors (the factor N is set by requiring canonical normal- 
ization), 



A. LEP Constraints 



una 




with A,B = {L, R} the chirality of the fermions. Typically A^ AB ~ 10 TeV, leading to mass bounds on 




(27) 
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FIG. 11: Tree level generation of multi-jet events at LEP due to flavor symmetric vectors. Similar diagrams (and 
bounds) exist for scalars and the first diagram is representative of the emissions of the vectors and scalars off the final 
state fermions. 



Some tree level processes that would lead to anomalous four jet events at LEP are shown in Fig. 11 These 
processes do proceed even in the absence of flavor violation. Furthermore, for vector pair production (the 
second diagram) the couplings of the Z or 7 to the vectors are fixed by gauge invariance. This implies a 
kinematic bound of My > 105 GeV on the flavor symmetric vectors. The analogous diagram for the scalar 
models is the origin of the kinematic bound on the scalar masses discussed in 0. 

The first diagram(s) can lead to a stronger bound for a vector or a scalar model as the vector or scalar 
is singly produced and has an order one branching ratio to anomalous multi-jet final states due to the flavor 
quantum numbers it carries. We find that this production process gives a cross section of 0.20|r/| 2 pb _1 for 
My = 120 GeV, or, with an integrated luminosity of 500 pb above this mass scale, a total of 100 1 77 1 2 events. 
The number of expected events falls rapidly with increasing My and is below 1 2 1 77 1 2 for My > 140 GeV. 
Assuming an order one branching ratio to anomalous multi-jet final states, this implies that the mass bound 
on the vectors from LEP is My > 150 GeV for 0(1) couplings. Colored vectors will have similar mass 
bounds from LEP from anomalous multijet events. The calculation for the scalar models is very similar, 
resulting in the same approximate bounds. 

In summary, we consider the minimum LEP lower bound for all vector and scalar models to be > 
150 GeV when order one couplings exist for the fields to quarks. A dedicated study would be required for 
each model to obtain more precise bounds. However, we note that while scalar (doublet) masses below 
150 GeV were only considered in model Si, see Eq. 22 for enhancing A l p B , this bound could easily be 



evaded by making the charged components of the scalar doublets heavier. Single exchange of the neutral 



components (the ones involved in tt production), as in the scalar analog of the first diagram of Fig. 1 1 would 
lead to kinematically suppressed virtual ti pair production. The LEP bound can be generically stronger for 
a color singlet that does not transform under the flavor group: for vectors I s — IV S the lower mass bound is 
> 1 TeV for order one couplings to leptons and quarks. 
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B. Electroweak precision tests 

We consider the models for massive Hp symmetric vector and scalar fields to be effective field theories. 
Up to this point we have assumed that the cut-off of the theory, A, is high enough that we can neglect 
operators suppressed by powers of A. For the scalar models, in principle, we do not have to consider effects 
of higher dimensional operators. In this case, for scalar models that are SUl(2) singlets, the results of Q 
show that the generic lower bound from fitting to EWPD is m s > 100 GeV. The LEP bounds we have just 
derived are stronger than these bounds for scalar models that do not break custodial symmetry. 

Below, we summarize our findings for model Si, for which particularly light scalar doublets were con- 
sidered in the context of Ap B . We find that it is easily consistent with S parameter bounds, even in the 
presence of nine light doublets, as in the Gf symmetric realization, rather than only two, as in the minimal 



Hp symmetric versions 1(a) - 1(c) (Eqs. ( |2Tj ) - ((23])) considered in the previous section (we have general- 
ized the analysis of ll37l for a single scalar doublet). The T parameter bounds are easily satisfied in the 
minimal versions, with neutral scalar doublet component masses as low as 110 GeV and heavier charged 
component masses of 150 GeV, a mass splitting that may be required in order to simultaneously evade the 
LEP bounds discussed above and enhance A l p B to observed levels. Mass splitting between the charged and 
neutral scalar doublet components can be obtained via Gf symmetric couplings to the SM Higgs vacuum 
expectation value. We have also checked the sizes of non-oblique corrections in the minimal versions of 
Si. In options 1(a) and 1(c), loop graphs containing the top quark lead to problematic shifts in the LH down 
and strange quark couplings to the Z. For instance, taking charged scalar masses of 150 GeV, neutral scalar 
masses of 110 GeV, and couplings r\ « 2/3 yields 5g SL /g SL , 5gd L / Qd L ~ —0.0046, roughly 2a in excess 
of the 2 a upper bound obtained, for equal shifts in the two couplings, from the Z hadronic width. In option 
1(b), loop graphs containing the charm and up quarks lead to a shift in the LH bottom quark coupling to the 
Z. For the same parameter choices (see Fig. [5]) the shift is 5gb L /gb L ~ —0.0028, which is consistent with 
Rb at 2a. We therefore conclude that whereas 1(b) is a viable option for enhancing Ap B , this is unlikely for 
cases 1(a) and 1(c). 

For the massive vector fields we do have to consider the effects of higher dimensional operators sup- 
pressed by the cut off scale as the theory is nonrenormalizable. This is particularly the case when con- 
sidering electroweak precision data (EWPD) for vectors. To the effective fields we have considered at the 
scale fi ~ My we should add dimension six operators. These arise from integrating out heavier modes 
of the UV complete theory. At the electroweak scale fi ~ mz we integrate out also the vectors and can 
match onto an effective field theory constructed from only the SM fields, appropriate for examining EWPD. 
The dimension six operators then receive contributions from vectors and from heavier fields from the UV 
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FIG. 12: Contribution to the operator 0\ generated by the SM through mixing with the SM B field, 
completion. 

For cases I s — IV S , the EFT will in general have a dangerous dimension 4 operator 0\ = V 1 (H* H). 
Integrating out the vectors gives the dimension six operator Ot = D ^ H\ 2 / My , well constrained 
by EWPD. For positive (negative) Wilson coefficient Ct the bound on the vector mass is My > 
5.6 (4.6) TeV x ^/\Or~\ (adding only this operator to the SM EWPD fit, and taking m h = 115 GeV) 11271 . 

The operator 0\ could be suppressed by the unknown dynamics in the UV theory. However, it is also 



generated via the quarks that the vectors couple to and SM interactions, see Fig. 12 Due to mixing with the 



SM B field (with coupling g) as shown an effective wilson coefficient of Ci ~ Sj y 2 rj 2 g 3 v A / (16tt 2 My) 2 
is generated, where yi is the quark Yukawa coupling and the sum is over all the allowed flavors (we neglect 
small Gf breaking). The bound is very weak as this is effectively a two loop effect due to mass mixing 
with the B field. This gives a bound My /rj 2 / 5 > 50 GeV for models II S , III S , IV S generated through SM 
interactions, weaker than the LEP bound discussed. The bound is even weaker if the vector only couples to 
down quarks and y^ is not enhanced from its SM value. Models that transform under flavor can still generate 
this operator through SM interactions but require the insertion of the appropriate number of Yukawas in 0\. 
If the vector transforms under SU(3) C the operator 0\ is forbidden. In general, if this operator is not present 
due to matching onto the underlying theory, the bound through custodial symmetry violating loop effects is 
weaker than the direct production LEP bound. 

Note that a key assumption of this analysis — fitting to EWPD with the SM supplemented with a 
single operator such as Ot — is not well justified. We leave a more complete analysis of EWPD to a 
future publication but note that non oblique observables, such as and the decay width of the SM Z 
boson potentially can place stronger bounds on the Hp symmetric fields with a large number of degrees of 
freedom. This is particularly the case when flavor off diagonal couplings are present and the top quark mass 
scale is present in large numbers of loop corrections to Rf, and Tz- 

C. Tevatron and LHC dijet Constraints 



At the LHC and the Tevatron the most relevant constraints on flavor universal couplings come from dijet 
resonance searches and dijet angular distribution studies. These constraints limit the couplings of the Hp 
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symmetric fields to the light quarks. In this section we give a brief account of these dijet bounds on scalar 
models Sy, Syi and vector models II , VI Q , and VI S , relegating details of the calculations to Appendix [B] 
These models are of interest for the tevatron's top quark forward-backward asymmetry, A l p B . The vector 
examples contain prominent dijet contributions in both the s and t-channels at the tevatron and LHC. Scalar 
model Sy is u-channel dominated, due to the flavor antisymmetry of its quark couplings (s-channel effects 
are present in the uc — > uc, uc — > uc channels); model Syi contains s-channel (uu — > uu) and -u-channel 
contributions. Related bounds on some of these models were also discussed recently in ll35ll36l[T28ll . 

The relevant scalar interaction Lagrangians in the Gf symmetric limit are given in Eq. ( |A37 1. The 
relevant vector Lagrangians are listed in Eqs. (|4])-(fT2|) in Section[IIJ and in Appendix [A] (with Gp breaking 
corrections also discussed). Expressions for the dijet cross sections in these models are given in Appendix 
|b| The new scalar or vector fields can have 0(1) couplings to first generation quarks, a necessary condition 
for obtaining large A t p B . Thus, they can mediate large contributions to 2 — > 2 quark scattering, potentially 
spoiling the good agreement between measurements and SM theory predictions for dijets. Tevatron and 
LHC bounds on their couplings to light quarks are collected in Table [V] They are derived using partonic 
cross sections evaluated at LO in a s , with the renormalization and PDF factorization scales identified with 
the transverse momenta of the two outgoing partons. Rapidity related experimental cuts on the two leading 
jets (the dijet measurements are inclusive) are emulated by imposing these cuts on the outgoing partons. 
Clearly, a parton level treatment introduces considerable uncertainties. Nevertheless, our bounds should 
still be indicative of the constraints one would obtain with more precise inclusive Monte Carlo simulations, 
carried out at NLO and including showering and full detector response. 

In Tables [V] and VI we quote upper bounds on the couplings -q and r/i (or f q ) of the light quarks to the 
scalar and vector fields, respectively, for a range of masses. There are four sets of bounds, obtained from 
searches for dijet resonances in the dijet invariant mass (Mjj) spectra at CDF HI 291 and LHC IU30H133l l. and 
from dijet angular distribution measurements at D0 [134] and CMS 1113511 . The CDF dijet resonance search 
M129II (corresponding to a luminosity of 1.3 fb _1 ) imposes a rapidity cut of \y\ < 1 on the two leading jets. 
For the LHC dijet resonance searches we use the CMS data [ 130] (2.9 pb -1 ) for 500 GeV scalar or vector 
masses, the ATLAS data IIT3TI (36 pb" 1 ) for 700 GeV masses, and the more recent ATLAS |[T32l (0.81 
fb _1 ) and CMS H 1 3311 (1 fb _1 ) data for heavier masses. These studies impose a rapidity cut of \y\ < 2.5 
(2.8 in H1321 ) on the two leading jets, with a rapidity separation for these jets satisfying | Ay\ < 1.3 (1.2 
in II132II ). Rapidity cuts are imposed in order to eliminate a large fraction of the dominantly i-channel, 
hence forward, QCD dijet background at large px', s-channel NP effects will be more isouopic. As already 
mentioned, we implement these cuts on the two outgoing partons. 

The dijet angular distribution measurements quote normalized differential cross sections 1 / Odijet da/dx, 
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Sy Mass 


TeV M j:j LHC 3/, , 


TeV X 


LHC X 


Syi Mass 


TeV Mjj 


LHC Mjj 


TeVx 


LHC X 


300 


1.0 




1.2 


1.1 


300 


0.3 


- 


0.4 


0.5 


500 


1.2 


n.b. 


0.5 


0.9 


500 


0.3 


2.2 


0.2 


0.5 


700 


2.0 


0.7 


0.7 


0.6 


700 


0.6 


0.2 


0.2 


0.2 


900 


2.5 


0.3 


0.6 


0.5 


900 


0.7 


0.1 


0.2 


0.2 


1100 


2.8 


0.4 


0.5 


0.6 


1100 


1.4 


0.1 


0.2 


0.1 


1300 


4.0 


0.5 


1.3 


0.6 


1300 


1.6 


0.1 


0.7 


0.1 


1500 


6.0 


0.6 


1.6 


0.3 


1500 


1.8 


0.1 


0.8 


0.1 


1700 


n.b. 


0.6 


1.8 


0.5 


1700 


2.0 


0.1 


0.8 


0.1 


1900 


n.b. 


0.6 


2.0 


0.4 


1900 


2.6 


0.1 


0.9 


0.1 


2100 


n.b. 


0.7 


2.1 


0.6 


2100 


3.0 


0.1 


1.0 


0.1 



TABLE V: Approximate upper bounds on the couplings of the scalars Sy, Syi to light quarks due to the measured 
dijet invariant mass spectra (labeled Mjj) and angular distributions (labeled x) at the tevatron and LHC, as explained 
in the text. The masses correspond to the scalar field flavors Sy and Syj =22=12 . If no bound is determined we denote 
this with "n.b.". 

where the angular variable is \ = (1 + 1 cos#|)/(l — | cos#|) and 9 is the scattering angle for the 2 — > 2 
partem scattering process in the parton CM frame. The differential cross sections are integrated over dijet 
mass intervals of a few hundred GeV in size, and results are presented in bins of x for < x < 16. 
The normalization, Odijct. is the measured cross section integrated over the dijet mass interval and over 
< x < 16- Whereas the x distributions for the QCD dijet background are relatively flat (due to t- 
channel dominance), more central NP effects, e.g., from s-channel resonances, will peak at low x- The 
D0 angular measurements II 1341 (0.7 fb -1 ) include a cut on yt>oost = 0.5 \yi + 2/2 1 < 1, where are the 
rapidities of the two leading jets. The CMS angular measurements [ 135 1 (36 pb _1 ) employ a similar cut of 
2/boost < 1.11- Again, we implement these cuts on the two outgoing partons. 

In general, the dijet invariant mass spectrum of a NP model can exhibit both an s-channel peak in the 
vicinity of the mass of the mediating field, and a monotonic rise relative to the SM prediction at larger 
Mjj. At the Tevatron, our scalar examples would possess the latter feature, due to u-channel effects, but 
an s-channel peak would be less prominent, particularly for model Sy- Therefore, for these models the 
bounds on the light quark couplings quoted in the Tevatron "Mjj" columns of Table [V] are obtained by 
requiring that the ratios of the predicted to SM da / dMjj spectra (with both evaluated at LO) lie within 
the PDF uncertainty band in Fig. lb of the CDF study 111291 . In particular, we do not make use of the 
bump hunter bounds on the new particle production cross sections. The widths of the intermediate scalar 
fields are identified with the two-body decay widths to quarks neglecting phase space, T = rf m s / (2 n) for 
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Ho 


VI 


VI S 


Mass 


TeV Mjj 


LHC Mjj 


TeV X LHC X 


TeV Mjj 


LHC Mjj 


TeV X 


LHC X 


TeV Mjj 


LHC Mjj 


TeVx 


LHC X 


300 


0.4 


_ 


0.9 


1.7 


0.6 


_ 


1.4 


2.2 


0.6 


_ 


1.7 


1.4 


500 


0.4 


0.8 


0.3 


1.3 


0.4 


0.9 


0.4 


1.5 


0.5 


1.0 


0.5 


1.4 


700 


0.4 


0.8 


0.3 


0.9 


0.4 


1.0 


0.5 


1.0 


0.5 


1.1 


0.6 


1.2 


900 


0.3 


0.3 


0.2 


0.7 


0.5 


0.3 


0.3 


0.9 


0.6 


0.3 


0.4 


1.0 


1100 


0.4 


0.3 


0.4 


0.8 


1.1 


0.4 


0.6 


1.0 


1.3 


0.5 


0.8 


1.2 


1300 


0.7 


0.5 


1.1 


0.9 


1.0 


0.6 


1.6 


1.2 


1.2 


0.6 


1.6 


1.2 


1500 


2.7 


0.5 


2.6 


0.9 


4.8 


0.6 


4.3 


1.1 


3.4 


0.7 


3.1 


1.2 


1700 


4.0 


0.5 


3.5 


1.2 


6.4 


0.7 


5.7 


1.6 


4.8 


0.7 


4.1 


1.6 


1900 


5.3 


0.7 


4.4 


1.0 


7.7 


0.9 


6.5 


1.3 


6.1 


0.8 


4.8 


1.4 


2100 


6.5 


0.8 


5.2 


1.4 


8.8 


1.0 


7.2 


1.8 


7.8 


0.8 


5.6 


1.8 



TABLE VI: Approximate upper bounds on the couplings f q of the vectors II Q , VI and VI S to light quarks, due to the 
measured dijet invariant mass spectra (labeled Mjj) and angular distributions (labeled x) at the Tevatron and LHC, as 
explained in the text. 



Sy and Syi- The vector examples can produce prominent s-channel peaks at the Tevatron. We therefore 



obtain the Tevatron Mjj bounds quoted in Table VI from the CDF bump hunter 95% CL upper limits on 



the product of a Z' production cross section x its branching ratio to dijets (Br) x acceptance (A), setting 
Br x A = 1 (see Table I of II129II ). For vector masses in excess of the last measured invariant mass bin, 
Mjj 6 [1225, 1350] GeV, we require that the ratios of the predicted to SM dijet cross sections in this bin 
lie below the 95% CL upper bound, obtained from the systematic and PDF uncertainty errors added in 
quadrature (the corresponding ratios in the lower bins would be much smaller than such bounds, given the 
widths of the vector fields). The intermediate vector two body decay widths to quarks are rj 2 my/(16 7r) 
for IIo and Vis, and rj 2 my / (96 it) for VIo. To obtain bounds on the scalar and vector couplings from the 
ATLAS and CMS dijet resonance searches we again use bump hunter 95% CL upper limits on new particle 
production cross sections, setting Br x A = 1. Specifically, for the CMS studies 111301 11331 we make use 
of the bounds provided for qq final states (see Table 1 of H130II for a 500 GeV resonance mass, and Table 1 
of H1331 ); for the ATLAS studies in II 13211 we use the upper limits in Table 3, and in B131H we use the bound 
in Figure 3 (for a 700 GeV resonance mass). 

To set the bounds on the light quark couplings from the dijet angular distribution measurements, we 
conservatively require that in each bin of x and in each Mjj interval, the addition of NP does not result in a 
shift of l/ffdijet dcr/dx that is larger than the la error. In the case of the D0 measurement 1113411 . the errors 
are statistics dominated (see the "EPAPS" document cited therein), and we can neglect the systematic and 
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SM theory errors. For the CMS study H135B we add the theory and experimental errors in quadrature. In all 
of our examples, the largest deviations occur at small and they tend to concentrate in the invariant mass 
intervals closest to the new particle masses, as would be expected for s-channel resonances. We can see 
from a comparison of all the bounds in Tables [V| and VI that the LHC bounds have superseded the Tevatron 
bounds for new scalar or vector masses > 1 TeV, and that they are particularly restrictive for model Syi- 

It is important to note that the constraints determined above are on the couplings of the Hp symmetric 
couplings of the new fields to light quarks. Only in the y t — > limit, when the full Gf symmetry is restored, 
are these couplings the same as the couplings involving heavy quarks. Thus, these bounds cannot be directly 
translated into limits on the allowed A FB , in general. Rather, these bounds dictate to what degree Gp — > Hp 
flavor breaking is required, in order to have couplings that are suitable for explaining the A l p B anomaly 
while being consistent with dijet constraints. Considering the benchmark scalar examples in Fig. [3j flavor 
symmetry breaking is probably not required in model Sy, while moderate breaking is required in model 
Syi- In the case of vector models VI Sj0 , we can see that in examples in which simultaneous agreement with 
all of the top data is obtained with the help of the Gf —> Hp symmetry breaking hierarchy f q f t << f q t, 
see Fig. [9] some flavor symmetry breaking would probably also be required in order accommodate the dijet 
bounds (f q < f qt ); however, in models in which agreement with the top data is achieved with the help of 
enhanced vector decay widths, see Fig. [2| little or no flavor symmetry breaking appears to be required by 
the dijet bounds. 

We close this Section with a brief discussion of collider constraints on scalar model Si. In particular, 
we consider the recent D0 bounds HI 1611 on anomalous dijet production in Wjj final states, prompted by a 
potential signal at CDF 11171 . and constraints from UA2 dijet measurements 111 18H . The scalar doublet fields 



Sij, i,j = {1, 2} couple to light quarks with equal Hp symmetric couplings rjij, see Eq. 18 which would 
lead to dijet peaks in Wjj final states via associated W + Sij production. D0 has presented 95% CL upper 
bounds on the cross section for such processes, ranging from 2.57 pb to 1.28 pb for dijet resonance masses 
of 1 10 - 170 GeV, i.e., in the preferred range for enhancing A l p B in model Si. We apply these bounds to the 
sum over cross sections for the four scalar doublet fields (charged and neutral components). The resulting 
95% CL upper bounds on the couplings rjij are w 0.2 for these masses. Finally, we find that the UA2 dijet 
measurements do not provide meaningful constraints on these couplings once finite width effects for the Sij 
and interference with the SM are taken into account. 
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D. Residual Constraints from FCNCs 

The fact that there is a global flavor symmetry Hp in the models we consider helps to avoid low energy 
constraints from Flavor Changing Neutral Currents (FCNCs) even though the NP has weak scale masses. In 
this section, as always when considering flavor bounds, we adopt a strict interpretation of the MFV hypoth- 
esis and assume that Gf is only broken by the SM Yukawas. The detailed analysis of FCNC constraints 
from B^s, D and K mixing is given in Appendix |C| while here we only collect the main results. 

We discuss the models in the two categories, universal and Yukawa suppressed, that we have introduced. 
(In the notation of Ref. |[T3ll these are the models that receive only contributions from class- 1 or class-2 
operators). The universal models give roughly the same contributions to and B s mixing amplitudes, 
when normalized to the SM. For vector or scalar masses about 1 TeV, the coupling constants have to be less 
than about 0.1. The models of the universal type are vector models III s o , IV Sj0 , VIII s o , IX s o and scalar 
model Sxiv- 

The Yukawa suppressed models give contributions to the mixing amplitudes that are CKM suppressed 
(with the same CKM structure as the leading short distance contributions in the SM). They are, in addition, 
Yukawa suppressed. The suppression can start at linear or quadratic order in light quark Yukawas. The 
FCNC bounds are satisfied for rj couplings of 0(1) for vector or scalar masses above ~ 1 TeV for linear 
Yukawa suppressed models. For quadratically Yukawa suppressed models the FCNC constraints are even 
weaker and are already satisfied for masses above ~ 300 GeV. 

Finally, meaningful constraints can be obtained for certain quark couplings in model Si from the B — > 
Kit branching ratios IT371I115I1 . For example, if the LH quarks in the Si Lagrangian of Eq. ( [T8] ) are defined 
in the up quark mass eigenstate basis, the couplings of the charged scalars in the quark mass eigenstate 
bases would include the terms 

Iki Ski ( V CKm) 3 kd' jL u' R + h.c, (28) 

where j, k = {1,2,3}. In the Gf symmetric limit, the contributions to B — > Kir, irir from tree-level 
exchanges of the scalars would exactly cancel. However, in the limit of large Gf — > Hp symmetry 
breaking (as would be required in order to maintain consistency with the Wjj bounds while accounting 
for enhanced A l p B ) the couplings would be constrained. The tightest bounds come from B — > K + ir~ . 
Borrowing from the analysis of OTTl . we find rjk l < 0.23 (ms/125 GeV), when considering the exchange 
of and separately. Thus, tree-level exchange of the scalar doublet fields could not account for 
A l p B in this case. However, if the couplings of the LH quarks in Eq. ( [T8j ) are defined in the LH down quark 
mass eigenstate basis, the scalars would only couple to bi in the Hp symmetric limit, as in Eq. 22 and 



40 



there would be no obstruction to A FB enhancement from B — > Kir. 



V. LHC SIGNALS 

We next turn to signals at the LHC from a flavor symmetric sector. The NP fields couple to light quark 
bilinears so that generically deviations in the phenomenology of dijet observables at LHC and the Tevatron 
are expected. No statistically significant deviation in these observables has been measured to date, thus 
placing bounds on the couplings of NP fields to light quark bilinears. For a set of sample models we derived 



the bounds in Section IV C showing that 0(1) couplings of these fields to quark bilinears are generally still 
allowed. Deviations from the QCD prediction can arise from an s-channel resonance effect, a t-channel 
driven rise in the dijet invariant mass spectra at large Mjj, or a combination of both signatures. 

A prototypical signal of colored Gf symmetric fields is pair production, g g — > V V and g g — > S S. 
This flavor diagonal scattering comes about due to the kinetic terms of the new scalar and vector fields, and 
the rate is driven by a coupling fixed by gauge invariance. The gg — > SS production cross section in the 
Gf symmetric limit as a function of scalar meson mass is shown in Fig. [13] As this signal requires the 
pair production of the NP states, there is a significant fall in the cross section with the mass of the scalars. 
A simple approximation to the fall in the pair production cross section as a function of ms is given by 
Ccff ~ (1, 16, 20) pb x exp(— m s /95 GeV) for the color representations (3, 6, 8). This approximates the 
full result within the PDF uncertainty band over the mass region 0.15 — 1.5 TeV. We expect the behavior 
of the g g — ^ VV cross sections to be of a similar order of magnitude. Note that the non renormalizable 
nature of the effective vector field will lead to corrections to the cross section that grow with ^fsjmy. These 
effects, however, are accompanied with the fall in the cross section with sfs due to PDF suppression. 

Dedicated four quark jet studies where one searches for pairs of equal mass resonances with no missing 
energy can place strong bounds on models of this form in principle. The cross section for 4-jets at the 7 
TeV LHC with all four jets having at least p T of 300 GeV (500 GeV) is 1.6 • 10~ 3 pb (1.6 • 10~ 4 pb), 
obtained using ALP GEN. Compared to the pair production cross section the QCD background is thus of 
moderate size. The condition that the four jets reconstruct resonances of approximately equal mass can aid 
in reducing the QCD four jet background. We leave a detailed study of the discovery prospects for a future 
study. 

There are a number of other LHC signal features that are common to Gf symmetric models. When 
a signal requires flavor violation, such as the production of same sign tops, the rate is suppressed by the 
required CKM suppression and Yukawas to generate the tops from the initial state - as in the SM. The 
i-channel models aiming to explain the A l BB anomaly without such a flavor symmetry can be strongly 
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gg -> SS (3,6) gg -> SS (8) 
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FIG. 13: Production cross sections of gg — > S S for color triplets (red left), sextets (blue left) and octet (grey right) 
scalar field S. We show a e f / = cr/ dim[-F] for a single flavor. We use MSTW208 pdfs and the PDF error bands shown 
are defined by taking fi = [2(2mg), (2ms)/2]. Details of the calculation are in the Appendix. 

constrained by bounds on same sign top signal at LHC H136B . The Gf symmetric models we have studied 
are essentially unbounded by this measurement. 

Other signatures that are linked to explanations of the A l p B anomaly do offer the prospect of constraining 
these models, such as a resonance decaying to t + j. This can be seen in ti + j production, and is a generic 
feature of i-channel models that can give enhanced A l p B not just Gf symmetric models. Prospects for early 
LHC discovery were explored in Ref. |[38ll44l . It is also feasible that for t-channel dominated models, the 
enhanced very forward ti cross section may be observable at LHCb HI 371 . We have shown that generically 
the new physics effect on the t i invariant mass spectra in the models we have considered are measurable in 
the near future. The cuts can be optimized at LHC to increase the prospects of measuring an excess of t i 
events due to new physics of this form while suppressing the SM background - see [ 138] for a recent study. 

VI. CONCLUSIONS 

In this paper, we have explored the collider physics of new physics sectors symmetric under the flavor 
group G F = U(3) Ur x U(3) Dr x U(3) Ql and its subgroup H F = U(2) Ur x U(2) Dr x U(2) Ql x U(l) 3 . 
These are the flavor symmetry groups of the SM quark sector in the absence of all yukawa couplings, and 
in the limit that only the top and bottom yukawa couplings are kept, respectively. Flavor symmetric fields 
cataloged in Tables I and II can have EW scale masses and 0(1) couplings to quark bi-linears. This is 
because the representations of these fields under the SM gauge group and the flavor group are motivated 
to address a generic experimental conflict - the flavor problem of NP. Moreover, the flavor symmetries 
eliminate NP contributions to single top and same sign top production, which otherwise would be tightly 
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constrained. For fields that transform under Gf and have 0(1) couplings to light quarks we find that a 
lower bound on their masses of roughly 150 GeV can be obtained from the absence of anomalous multi-jet 
events at LEP2. More detailed studies of EWPD, however, have the potential to improve these constraints. 

Upper bounds on the couplings of Gf symmetric fields to light quarks can be set through studies of 
dijet resonance searches and angular distributions at the Tevatron and LHC. Based on the partonic LO study 
of dijet production we have determined these bounds in a number of sample models. Our results indicate 
that the LHC data sets are starting to constrain the allowed couplings to be < 0.1 in some of these models, 
in particular if they lead to resonant uu quark pair production as in models with color sextet scalars (e.g., 
model Syi in Table [II]). 

FCNC bounds on flavor symmetric models depend on how Gf is broken. In the Hp symmetric limit 
there are no NP contributions to meson mixing. We have analyzed NP contributions to meson mixing 
assuming MFV breaking of the Hp symmetry. The contributions to meson mixing are of second order in 
the product of MFV flavor breaking coefficients and vector or scalar couplings to quarks. We find that 
for models that contribute to Bj and B s mixing with the same product of CKM elements as the SM box 
diagrams (and no additional light quark yukawa coupling factors) these coefficients have to be < 0.1 for 
vector or scalar masses of 1 TeV (these are vector models III s o , IV S)0 , VIII s o , IX s o and the scalar model 
Sxiv)- The bounds on couplings for the remaining models are much looser. 

The fields we have introduced and studied are of particular interest in the LHC discovery era. As 
an application we have examined the potential impact of such field content on two current experimental 
anomalies, the top quark forward-backward asymmetry, Ap B , and the B St d dimuon charge asymmetry. 
Many of the flavor symmetric models can yield large A t p B , while maintaining agreement with the total 
top pair production cross section. However, simultaneous agreement with the m t i differential cross section 
is more difficult to achieve. Generically, scalar models have difficulties accommodating all of the top 
quark data. The scalar models can either give u-channel contributions to tt production (models Sy — Sxiv) 
or t-channel contributions (models Si — Siv)- The u-channel models exhibit significant tension with the 
differential spectrum. Of the t-channel models we have checked that model Si shows good agreement with 
the differential spectrum, for masses of order 100 GeV. While it has difficulty reproducing the CDF m t i 
dependent Afb at the one sigma level, it should be noted that the recent D0 A FB measurement appears to 
exhibit a smaller m t j dependence. 

Vector models can contribute to ti production in the u-channel, s-channel, t-channel, or the s-channel 
and t-channel simultaneously. Phenomenologically preferred t-channel dominance can be a feature of the 
flavor structure (models VH S)0 ), or can be attained either through some flavor breaking or additional decay 
channels to broaden the s-channel decay widths. This leads to good agreement with all the top quark data, 
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with a preference for lighter vector meson masses (see models VI Sj0 in Fig. [9j or Fig. [2]and Ref. 051 . and 
models VII Si0 in Fig. [8] ). Note that flavor breaking or broadening of the widths (via non-qq final states) 
also insures consistency with resonance searches in the dijet data. 

When studying the potential of new physics explanations of the A l p B anomaly we have determined 
model independent acceptance corrections. These are appropriate for comparing against Tevatron data that 
can be applied to any model to good approximation at the partonic level- these results are given in Section 
MIA II 

Assuming MFV for Hp symmetry breaking, we have systematically examined the effects of scalar and 
vector fields on the dimuon charge asymmetry anomaly. Flavor symmetric models that can easily give 
relatively larger contributions to B s mixing than to mixing are preferred. Several flavor symmetric 
models can accommodate this pattern with tree level exchanges: vector models I s o , V s>0 and scalar models 
Srv, Sviii, with 0(1) MFV flavor breaking coefficients and NP masses of no more than a few 100 GeV, and 
models Xg 6 , Sin, for masses around 1 TeV. None of these, however, can also explain the large A l p B , while 
simultaneously obeying all constraints. The dimuon charge asymmetry anomaly can also be explained by 
the experimentally less favored possibility - equal contributions to B s and Bj mixing amplitudes (once 
normalized to the SM contributions). There are a number of models that can contribute to ti production as 
well as B Sj d meson mixing. Of these, vector models VII so and scalar model Si are preferred, leading to 
both very small deviations in the differential ti spectrum and enhancement of Ap B . The contributions to 
B s d arise at 1-loop and can be of the right size. 

The models we have introduced and motivated can be discovered at the LHC through a number of 
processes: ti + j with t + j forming a vector or scalar resonance peak; deviations in the m t j differential 
cross section; pair production via gg — > SS, VV, where two resonance peaks are reconstructed in four jet 
final state; and deviations in dijet observables. 

In conclusion, flavor symmetric fields by virtue of their flavor safe properties can have masses near the 
electroweak scale. They are therefore candidates for discovery at the LHC, as the lowest lying states of 
more complete models, which may or may not be directly linked to electroweak symmetry breaking. Flavor 
symmetric sectors have a rich, predictive, and mostly unexplored phenomenology. It is interesting to note 
that the list of fields that are Gf symmetric and can directly couple to quark bilinears is relatively short. 
Systematically exploring the phenomenology of flavor symmetric sectors is therefore feasible. 
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Appendix A: Flavor Symmetric Lagrangians 

In this appendix we collect the Lagrangians of the vector fields in nontrivial Gf representations. The 
form of the scalar Lagrangians is given in |f3l and we use that notation in the main part of the paper and 
in the Appendix for ease of comparison (for concreteness we also show explicit form of the relevant scalar 
Lagrangians at the end of this Appendix). For parameterizing the breaking of Gf we use the language of 
MFV, where the breaking is due to SM Yukawas only. For the unitary transformations that diagonalize the 
quark mass matrices we use the following notation 

U(u,R?YuU(u,L) = V2M u /v = Y^ g =± diag(0, 0, y t ), 



(Al) 

~. .-!;.,.-. ;n ii ,,.\ 

These relate the flavor (qi) and mass eigenstates (q0 



U(d,R)^Y D U(d,L) = V2M d /v = y^ iag ~ diag(0,0,y b ) 



u L =K(u,L)u' L , u R =U(u,R)u' R , 
d L =U{d, L)d' L , d R = U(d, R)d' R . 



In ( |A1| ) M. u and Aid are the diagonal up-type and down-type quark mass matrices and v ~ 246 GeV is the 
SM Higgs vev. The Cabbibo-Kobayashi-Maskawa (CKM) matrix Vckm is then given by 

^ckm =U(u,L)^U(d,L). (A2) 

Two basis choices that we will frequently use are the "up" and "down" quark mass bases. In the "up" 
quark basis Y v -»• (Yjf ag ) 2 , Y ] D Y D -> V C km(Y^ s ) 2 V^ km , so that down-quarks are rotated from the 
mass eigenstates to keep interactions with the W bosons diagonal. Similarly in the down quark mass basis, 
the up quarks are rotated by V ] CKM , so that y} Y v -> ^(C^^ckm, ^ Y D -> {Y^ f. 



1. Vector fields I - IV 



The vector fields I s o — IV SiQ are flavor singlets, but can be either singlets or in adjoint representations 
of color and/or weak SU(2)l. The interaction Lagrangian with quarks then has the same forms as the 
interactions with corresponding gauge fields would have. For instance, for I s o , it is 

d = r,d R fd R , (A3) 

where = V^T A for I Q fields, with T A the Gell-Mann matrices normalized to Tr(T A T B ) = 5 AB /2. 
The interaction Lagrangians for II S]0 , III Sj0 are obtained from the above through replacements, dR — > ur 
and dR — > Ql, respectively. In the interaction Lagrangians IV Sj0 the short-hand notation for the sum of 
vector field multiplets is = V^a a /2 for IV S model and = V^' a T A a a /2 for IV G model, with a a the 
Pauli matrices. 
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2. Fields V s 



These fields are color singlets (octets) and SU(3)d octets. The interaction and mass terms in the La- 
grangian are similar to the case of the SU(3)u octet fields VI Sj0 given in Eqs. @-(|2]> and Eqs. ( |A4[ )-( |A6] ), 



but with the obvious u R — > d R and y t — >• yb, Yjj — > Yd replacements. In the SM, y^ ~ 0.02 so that 
all the vectors in the multiplet are degenerate to a good approximation. The degeneracy can be lifted in 
theories with more than one Higgs where y^ can be 0(1), as occurs for example in the large tan (3 limit of 
the MSSM. 



3. Fields VL 



These real fields are 1, 8 under color and 8 under flavor. We denote them by V B and V^' B , with color 
(flavor) labels A(B) = 1, . . . , 8. The quark-vector interaction Lagrangian and the mass terms were already 
given in Sectionjn using the shorthand notation V* = T B V B , V° = T A T B V^' B . Here we write out the 
terms keeping all the indices explicit. The Gf symmetric quark- vector interactions Q are 



Aa s = Vl (urT* (T B M B (u R ) atj , £ V i = (0*)°'' (T A f a {T B )if A > B \u R ) PJ 



A \P (T B VW A > B t 



(A4) 



where a,/3, 7... (i, j, k...) are color (flavor) indices. The couplings rjl'° are real because the interaction terms 
are hermitian. The Gf — > Hp breaking interaction terms due to Y v Y^ insertions ([5]) have the following 
explicit form in the mass eigenstate basis (keeping only the terms proportional to yt) 



A£vi = yt (V2nu> R rf A > B (T A f a {T B )\{u' R )^ + h.c. 



^(n' R r% A ' 8 (T A fM,, (A5) 



AC 



vi s = Vt 



{u R rf B {T B Uu' R U + h.c. 



A 



y ^(u' R rr* (T*)Uu' R ) a . 



(A6) 



4. Fields VIL 



These fields are in the (3,3,1) representation of flavor. As for fields V, VI we use the short-hand notation 



V s = T B V B 



xro _ q-ArpByA,B 
v fi ' v fj, ) 



(A7) 



with T A the color Gell-Mann matrices. The flavor matrices T B now also include an identity matrix, so that 
for B = 1, . . . , 8, T B are the Gell-Mann matrices, and T 9 = 1/^6. Note that vjl / V^, because the 
fields are in the (3, 3, 1) representation of Gf- The transformation to the mass eigenstate basis is 



V? = U{d,R)V"U{u,R)\ a = o,s, 



(A8) 
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explicitly showing that V^' s are not hermitian (note also that V B , V A ' B carry nonzero hyper-charge). The 
tree level quark coupling Lagrangian terms are (no summation over a = o,s) 



£vn a = rfid R YZu R + h.c, 



with rft s complex in general. Explicitly the two terms are 



(A9) 



£vn = ^(V/' B T B )j(4) aj 7"(r A )fW ftll +h.c, 
£vn s = r, s 1 (V B T B y j (d R )^^(u R ) a , i + h.c.. 



(A10) 
(All) 

The flavor breaking introduces corrections to these interaction terms from insertions of the Y (J and 
Y D Yp. Keeping only the terms proportional to y 2 for Gf — > Hp, the leading contribution is given by 
V^dRYuYjjYyiiii, which gives in the mass eigenstate basis the following interactions 



AC 



VIIo 



Vt 



V2(d' R r i 7 M (T A t $ R )p, (V, A > B T B )f + h.c, 



A£ V ii s = V 2 t [v2(d' R )^^(t' R ) a ] (V B T B )f + h.c. 



(A12) 
(A13) 



The mass terms are 



£vn a = 2(1 + 5 a>0 ) [m 2 Va Tr [vj V^] + A (#t H) Tr [v£ V^] } , 



a = o.s, 



(A14) 



while the top quark Yukawa splits the mass spectrum through corrections to the mass Lagrangian of the 
form 2(1 + 6 a:0 ) CiTr[V^Y v Y^ V^]. The mass of the field V^ 8 gets split from the rest by the term 2(1 + 
<5<i,o) C2Tr[V^Y u Y^ [ V^ ajf Y D Y^]. These splittings result in the following mass spectrum 



2 2 2 2 i i , 

m 1 = m 2 = m 3 = m v I 1 + 



m A = m 5 = m 6 = m 7 = m v [ 1 + - + y Vt 



2m 2 v J ' 
A v 2 + Cj 



m| = my I 1 + 



A V 2 2Cl 2 ^ 2C 2 2 2 

2^ + 1^ + X^ 



A ^ Cl..2 , C2 _2 „.2 

3 



m 9 = m v [ 1 + -— + + ^ yl y{ 



(A15) 
(A16) 
(A17) 
(A18) 



where we neglect terms suppressed by off-diagonal CKM elements or first two generation Yukawas. 



5. Fields VIIL 



The scalar matrix of fields V£ = T B V B , and the octet matrix of fields V° = T A T B V A ' B , now trans- 
form as a bi-fundamental of U(3)q, i.e. V^' s — > VqV/?' s Vq for a transformation that takes the left-handed 
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quark fields to Ql — > YqQl (suppressing the flavor indices). Note that the fields V^ ,s are hermitian, 



v; 



o,s\ 



V^' s . The interaction with quarks is then given by 



£vm a = vIQlYZQl 



(A19) 



or writing out the color and flavor indices explicitly 



Aan = nlV^{Q L y> % l»{T A y a {T»)l{Q L ) M , 



tvm. = v{V*{QLT> i r(T B ){{Q L ) 



(A20) 
(A21) 



The flavor symmetric mass Lagrangian for the vectors is 



/"mass 
^-VIIIo,.. 



(1 + S a>0 ) [m 2 v Tr 1>] + A (ift H) Tr [>„ } 



(A22) 



where the color and flavor indices are suppressed. The masses are split by the presence of the SM Yukawas 
breaking the flavor group. In contrast to cases V and VII here both insertions of Yfi Y v and Y D are 
possible. The flavor breaking terms are 



V,YlY D V» 



+ C 2 Tr 



Y' Y V Y' Y 



+ ---(A23) 



with all real (complex couplings are possible at higher orders in Yukawa insertions). In general, the 
masses for are not diagonalized by the U(u,L) or U(d,L) unitary transformation, but by a unitary 
transformation that differs by ~ Vckm from the two. For simplicity, we work in the limit where the con- 
tributions from yi, are much smaller than from y t . In this limit the vector field mass matrix is diagonalized 
using U(u, L), giving a spectrum 



2 2 2 

m l = m 2 = m 3 



2 2 2 2 

m 4 = m 5 = m 6 = m 7 



TUy 



TUy 



rn 



V 



1 + 



1 + 



1 + 



2 m\ 
A v 2 
2 m v 

Xv 2 

2^ 4 T"T+.ru, 



2 

Cl 2 , C2 ,\ 

— yt H — "' 

3 yt 3 



(A24) 
(A25) 
(A26) 



The mass degeneracy is lifted at y 2 order from a contribution of the form y 2 (Vckm^C^km^ (i.e., from 



the term proportional to ([ in ( A23 1). 

We are now ready to write down the flavor breaking couplings between the vector and the quark fields. 
These are given by 



A£vm a = [v^QlY^Y^Ql + h.c] + [^Q L Y ] D Y D f^Q L + h.c] 
+ [vtQiYiY^Y^QL + h.c] + • • • , 



(A27) 
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where we only display one of the possible terms with four insertions of Yukawa matrices. Note that 772,3,4 
can all be complex. We first display the couplings to the up-quarks in the "up" basis 



A-Cviil, 



( u L)il3 



(h) a {T A f a {T u )\% 

+ yfal{V CKU ?i (« [^L) ia (Ty a {T B ))Y^ B (u' L ) kp 

+ Vhhl (^CKM)f (^CKM) J 3 [(?LTtT A t {T B )W*' B (u' L ) 3 p 
The couplings to the bottom quarks we display in the "down" basis 

A£ VIIIo = yhmKufi (^ckm)! \{d' L r(T A t(T B )]Y^ B (d' L ) 



y L ) a {T A t{T B )\ rf> B (d' L ) 



)k/3 



(A28) 



(A29) 



+ yhfa2(vS KM )nvcKM)l [(d' L ) ia (T A t(T B )*r A > B (b' L )p . 

For each of these two cases, the vector fields have to be put in the "up" or "down" basis respectively. In the 
limit where we neglect insertions, the "up" basis coincides with the mass basis for the vector fields, while 
the vector fields are not mass diagonal in the "down" basis. The couplings for the singlet case are obtained 
from the above expressions by replacing (T A )a — > $a- 



6. Fields IX s o 

The discussion for this case is very similar to the previous one. The fields are given by V B ' 1 for the color 
singlet and V A ' B ' 1 for the color octet case, where i = 1, 2, 3 is the weak SU(2)l index. All basis indepen- 
dent expressions apply also in this case, if the matrices of fields are defined as V* = T B a i V B ' i /2, and the 
octet matrix of fields V° = T A T B '^V^' 3 ' 1 / '4, where o % are the Pauli matrices. Phenomenologically the 
most important difference is that there are now charged currents. 



7. Fields X 3 6 

These fields are weak doublets and belong to the bi-fundamental representation of the flavor group 
(1, 3, 3). They are either anti-triplets of color, in which case we have the fields (V^ t )^, or they are sextets 
of color, in which case the fields are {V^- a p = (V^)l • a a . Here r is the weak SU(2)l index, a, j3 are 
the flavor indices, while i and j are the indices of the (1, 3, 1) and (1, 1, 3) representations respectively. In 
Section|IT]we have written the Lagrangian in short-hand notation, while here we show the terms keeping the 
indices explicit. The tree level quark coupling Lagrangian terms are (suppressing the SU(2)l indices) 

C Xs = Vi e«/3 7 (d R ) ai 7 M {Qlf j + h.c, (A30) 

Cx 6 = Vi (yjij^p (d R ) ai ^ (Ql)^ + h.c, (A31) 
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which in the mass eigenstate basis are 

Cxi = m ^ {V^J. (d' R r l 11 {u'ifi + m e aM (Vckm)* (V?)J fc K*) Q * V (d'l)^ + MA32) 

c X6 = ??i(^%>,M^r^^ (A33) 

with (V^ 1 ) and (V^) the charge and neutral vector fields, respectively. 

8. Fields Xlg 6 

The analysis is very similar to case X. The couplings to quarks, vector fields mass spectra and splittings 
can be obtained from the discussion on case X with the replacements d R <h> u r , Yd <h> Yjj. For instance, 
the quark coupling Lagrangian terms are (suppressing the SU(2)l index) 

Cx h = m ^ 7 WJ- (u R r ^ (QlY 3 + h.c, (A34) 
£xi 6 = m (^)m>,/3 {u R ) m ^ (Qlf 3 + h.c. (A35) 

9. Scalar Fields I, V, VI, IX, X 

For convenience, the interaction Lagrangians for the scalar models under discussion, in the Gf symmet- 
ric limit, are provided below: 

Ci = r] ( Si j u iL u jR + S'jdiLUjn) + h.c , (A36) 
C y = T,<?fre ijk Sj x v? Rfj u k Ri + h.c, £ yi = V S a > b (f a Ut b r^R a <p + h.c, (A37) 

C lx = r ] e a ^Sr j d i R0 u i Rl + h.c, C x = r, S^f^ d Ra u^ + h.c, (A38) 

where the rj are dimensionless flavor universal couplings, Latin (Greek) indices label flavor (color), and 
Lorentz spinor indices are suppressed. The generators T a and T h of the symmetric sextet flavor and color 
representations are the symmetric 3x3 matrices, normalized such that Ti(T a T b ) = 25 ab and a,b = 
1, • ■ ■ ,6. The scalar fields appearing in £i,v,vi,ix,x are canonically normalized. In C\ we show explicitly 
the Yukawa couplings of the charged and neutral components of the scalar doublets. The relevant NP cross 
section formulae for top quark pair production can be found in Appendix B. The decay widths of the scalars 
are nrj 2 ms/167r, where k = 1, 4, 2, 1, 1 in models Sy, Syi, Six, Sx, Si, respectively (assuming all quark 
decay channels are open and ignoring phase space effects). 
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Appendix B: Details on 2 — > 2, 3 scattering calculations and phenomenology 

In this Appendix we give some details of the calculations for tt and light quark qq production in the 
scalar and vector MFV models. The comparisons with experimental data at the Tevatron and LHC were 
done in Sections UTTland llV CI 



a. 2—^2 scattering in the scalar models 



The relevant interaction Lagrangians for the scalar models are given in Eqs. ( |A37| ) - ( |A36| ). The general 



formula for the spin averaged matrix element squared for qq — > tt, including interference with the SM, for 
the exchange of a color triplet or sextet in the n-channel is |[T8l 

~9sV 2 C {{sm 2 + u 2 )(u-m 2 s )\ 4f? 4 C 2 u 2 



.s-^lqiv Cn ( is mi + Uf)(u - m 2 )\ 

M 2 =My ^ -(V-* 9X9 T.9 9 1 + 

1 4^ [ s \ (u- m 2 ) 2 +T 2 m 2 J 



,2 



(Bl) 



(u — m 2 ) 2 + T 2 m 2 

where vt = v — m 2 for u = {s, t, u}. For the exchange of a color singlet or octet in the t-channel one has to 
make the replacement u—>t above. The color factors for an exchange of a (1, 3, 6, 8) of SU(3) C are given 
by C = (4, 1, -1, -2/3) and C 2 = (9, 3/4, 3/2, 2). 

The differential partonic 2 — > 2 cross section averaged over initial spins and colors and summed over 
final states is given by 

da 



dz 32 tt s 



Y^lMij] 2 . (B2) 



where £\ . |.Mij| 2 = \ \ Ylij \-M-ij\ 2 . Here (3 = y/l — \m 2 j s is the velocity of the final state quark (with 
mass m) in the initial state partonic rest frame while z = cos 9 and 9 is the scattering angle in the partonic 



CM frame. The sum in Eq. (B2i is over all contributing matrix elements. The initial states are weighted 



with the appropriate PDFs to obtain the hadronic cross section. When evaluating the dijet constraints we 



use MSTW2008 PDF sets [94]. The renormalization scale is set to the p T of the jets, fj, = Mjj \[\ - z 2 /2, 
where z = cos 9 is the partonic scattering angle from the beam line and Mjj is the invariant mass of the 
dijets. 

In the case of light quark pair production, in models Sy and Syi there are n-channel exchange con- 
tributions to u u — > cc,cc — > uu, and in Syi there are additional n-channel exchange contributions to 
uu — > uu, cc — )• cc. The corresponding matrix elements squared can be obtained from Eq. ( |B1[ ), taking 
care to include factors of 2 and 4, respectively, in the first and second terms for the latter subprocesses. 
s-channel exchange in models Sy and Syi contribute to uc — )• uc + cc, and in Syi there are additional 
contributions to uu — > uu + cc, cc ^ cc + cc. The corresponding cross sections can be obtained via 
crossing symmetry from the above, and including t channel SM contributions, also see |[36l . 
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Following M139II , experimental cuts on the rapidity y and the rapidity separation Ay of the two leading 
jets in dijet production are implemented at the partonic level in the 2 — > 2 subprocesses via 



da 

dM. 



jj 



on^I dy B V/ i (T e ^)/ i (Te- 



Vb 



da 

dz 



(B3) 



where ajsr = 0.3894 • 10 9 pb, s is the hadronic center of mass energy, r = y Mjj/s, fij are the MSTW 
parton densities, and Yb is the boost rapidity of the subprocess frame. The rapidity cuts 1 2/1,2 1 < ^max 
correspond to the integration limits 



y° B = Min[y max , - log[Af£./s]/2] 



(B4) 



and a cut Ay < (A?/) max enters the angular integration limits as 

z° = Min[z max ,tanh(y max -| yi j| 
where z max = (Exp[(Ay) max ] - l)/(Exp[(Ay) max ] + 1). 



(B5) 



b. Vector models 2 — >■ 2 production 

The vector models can be divided into two sets: the cases I-IX that couple only right-handed or left- 
handed quarks, and the models X and XI, where the right-handed quarks are coupled to left-handed (charge 
conjugated) quarks. To obtain a general equation for the qq — > tt cross section for the first subset of models 
I-IX, we consider a generic Lagrangian (with trivial changes in notation if vectors couple to right-handed 
quarks) 



£ = MlYql + fthftL + [f q tq L ft L + h.c], 



(B6) 



where we have suppressed flavor, weak and color indices. This gives for the weighted average of the 
amplitude squared for the qq — » ti scattering 



\M\< 



V" V V 



Amy 



+ 



\U? f q ft(s v tv + r 2 v m 



18 (s v + r 2 v m 2 v )(t 2 v + r 2 v ml 



s 2 [ (l + /3 e ) 2 + 2 



ml 



sm$ 



while the interference with the LO SM one gluon exchange gives 



SM) 



+ 



anas 
9 

87ras 
9 



C 4 |/, 



t V 



t 2 v + Tim 



V""V 



r fgftsy 

C5 4 + r? / m? s 



V"°V 



s 4 



16- 



m 



(B7) 



(B8) 
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Case 


Couples to q 




h 




ft 




fqt 


C{ 




ns ns ps 
u 3 u 4 u 5 


ci 


no 
°2 


no 
°3 


no 


no 


Is.o 


d, s, b 




m 












9 














2 











2 


Hs,o 


u, c 




m 




vi + vim 







9 














2 











2 


III S , 


u, d, s, c, b 




m 




m + VtV2 







9 














2 











2 


iv 8 , 


u, c 




m 




Vi + VtV2 







9 


-6 


3G 


8 





2 


4 
3 


8 


4 
3 


2 




d, s, b 








m + vim 




9 


-6 


3G 


8 





2 


4 
3 


8 


4 
3 


-2 


V s , 


d, s, b 




m 












i 

4 


i 

4 


9 

4 


2 





l 

is 


1 

18 


l 

2 


1 

3 


l 

3 


VI S , 


u, c 




m 




+ 2y 2 Re{ m ) 




+ vim 


1 

4 


1 

4 


9 
4 


2 





l 

18 


1 

18 


1 

2 


1 

3 


1 
3 


VII S , 


d, s, b 












Vi 


+ vim 








9 
4 


2 











1 

2 


1 

3 





VIII S , 


u, d, s, c, b 


m - 


V y 2 2Re{r, 2 )5 q . b 


m 


+ y 2 2Re{ m ) 


m 


+ vim 


1 

4 


-1 

4 


9 
4 


2 





1 

is 


1 

18 


1 
2 


1 

3 


1 
3 


IX S:0 


u, c 


m - 


Vy 2 2Re{r, 2 )5 q , b 




+ y 2 2Re{ m ) 


m 


+ vim 


1 
4 


1 
4 


9 
4 


2 





l 

IS 


1 

IS 


1 

2 


1 

3 


1 
3 




d, s 


m - 


V y%2Re(m)S q , b 


m 


+ y 2 2Re{ m ) 




+ vim 


1 

4 


1 

2 


9 


4 





1 

18 


1 

9 


2 


2 
3 


1 

3 




b 


m - 


V y 2 2Re{j] 2 )5 q . b 




+ y 2 2Re{ m ) 


m 


+ vim 


1 


2 
"3 


4 


8 
3 





2 
9 


4 

27 


8 
9 


4 
9 


2 
3 



TABLE VII: The contributions to ti cross section from MFV vector resonances. The coefficients C,; are to be used 
in Eqs. jB7\ , ( |B8| ), where we have kept for simplicity only leading flavor breaking in the couplings of quarks to the 
vector, while neglecting the flavor breaking in the vector propagators (and in models IV, IX, X, XI also neglecting the 
electroweak symmetry breaking contributions to vector masses). 



Here sy = s — my, and similarly ty = t — m v ,uy = u — m v . Explicitly one has 



t = rrif s(l 



u = mf - ^s(l + /3g), 



(B9) 



with f3g = f3 cos 6, while j3 = y/l — Am^/ s is the top velocity in the it rest frame. We have used the same 
notational conventions as in the scalar case above. This formula should be summed over all possible initial 
quark states q, which depends on the particular MFV vector field considered. 

The expressions for the contributions to the ti cross section for the second subset, cases X and XI, are 
very similar to the above. The important difference is that both left-handed and right-handed quarks couple 
to the vectors at the same time. The general form of the interaction Lagrangian in this case is 



C = I^UlPl + f R P R )Va q c + h.c, 



(BIO) 



which gives for the weighted average of the amplitude squared contribution to qq — > tt scattering 

12s 2 



Y,\M\ 2 = C[(\f L \* + \f R \ 



36lly 



^((l + /3,) 2 + 16^)+(l-/3 e ) 2 



+ 



+ C' 2 2\f L \ 2 \f R \ 2 



12s 2 
36u v 



in 



2, rc v) +41 _ 2 t 



. m 
s 



2\ 2- 



(Bll) 
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The interference with the SM one gluon exchange is 

Y2Re{MM* SM ) = ^7ra s C 3 (\f L \ 2 + \f R \ 2 ) — 

* — ' y Uy 



11 4 1 



(1-^ + 8^1 + ^(1 + ^)2 + 4."^ 



(B12) 

For the coefficients in these equations we have C[ = C' 2 = 1(1/2), C' 3 = 2(1) for 3(6) cases, while 
fL = m + Vt V2> fn = for model X and f L = r)\ + y 2 7] 2 , /r = m for model XI. 

For light quark pair production we use the generic form of the interaction Lagrangian ( |B6| ), but replace 
the top quarks with light quarks. We denote by f q the flavor diagonal coupling of vectors to quarks (these 
will contribute in the s-channel). The flavor off-diagonal couplings for terms contributing in the i-channel 
are denoted as f q q > . In the numerical analysis we set f q = f qq i , which follows from Hp symmetry. 

We examine dijet constraints on three sample models, II Q , VI S and VI Q . The spin and color averaged 
amplitudes for uu — > cc or cc — > uu processes are 

^\u\* = \Ml Ciu2 i CzfMsytv + TWv) 2, IAc| 4 C 3 u 2 



^ ' 9 s 2 v + T v m 2 v 9 (s 2 v + Tlm 2 v )(t 2 v + Tlrn 2 v ) 9 t 2 v + T 2 v m 2 v 

(B13) 

and 

^2Re(MM* SM ) = (C,\U 2 % 2 + C 5 /"^ 2 ) ~, 014) 

where i = — s(l — cos#)/2, u = — s(l + cos6>)/2. The Cj's are 



for model VI Q , 



for model VI S , and 



Ci — — , Co — — . C 3 — — , Ca — — , Cr — — (B15) 
18' 18' 2' 3' 3 



Ci = i, C 2 = -i C 3 = ^, C 4 = 2, C 5 = 0, (B16) 



Ci = 2, C 2 = 0, C 3 = 0, C 4 = 0, C 5 = 2 (B17) 



for model II . 

For the processes uu — > uu or cc — )■ cc the spin and color averaged amplitudes squared are 

Y\M\ 2 = f 4u 2 ( „ % _ + 2C 2 , _ t V 2 tV 2 + Jt m L 2 , + 2 ^ T ) (B18) 
^ 9 V4 + r y m y (s 2 v + T 2 v m 2 v )(t 2 v + T 2 v m 2 v ) t 2 v + T 2 v m 2 v J 

We have to add new terms to the interference between the NP and SM, because now there is also a i-channel 

contribution in the SM. The modified equation is 



^2Re(MM SM ) - — u f q K - [t + T ) + s 2 + f^m 2 { V + T 



(B19) 
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The d 's are 



p p ps pt ^ ps ^ pt 

1 — g) 2 _ ~27' 3 ~~ g' 4 ~~ _ n' L 4 — qJ <-5 " ~- l ~ 



2 

3' 



;> "-5 



in model VI Q , 



1 

3' 



C\ — 1, C 2 — -, C3 — 1, CJ = -, C4 — 0, C5 — 0, C5 



4 
3" 



in model VI S , and 



2 

3' 



/"7 — p r) pS — pt o pS r) /it 

1 — Z, L 2 — — -, C-3 — Z, C 4 — — — , C 4 — Z, L 5 — Z, L 5 



2 

3 : 



2 
3 



(B20) 



(B21) 



(B22) 



for model II G 

Using crossing symmetry, the processes uc — > uc + c.c. are obtained from Eqs. ( |B13| >, ( |B14| ) via the 



substitutions s — > t, t — > s, u — » u; the processes uc — ^ uc + c.c. are obtained from Eqs. ( |B13[ ), ( |B14| ) via 
the substitutions s — > t, t — > u, u — > s; and the processes uu ^ uu + c.c, cc ^ cc + c.c. are obtained 
from Eqs. (|B18|>, (|B19|) via the substitutions s — > t, t — > u, u — > s. 



c. LEP e+ e 



f fV production 



Consider a vector Lagrangian of the form C = rj ipR 7^ i/jr. Then the amplitude squared for the 
production of a color singlet vector V through e + e~ — > 7* — > Vf f is given by 



\M\< 



8e 4 Q'N c \ V \ 



2 a 2 



StS 



1°2 



-2 (2ml(s 1 s 2 - 3xti - s 2 t 2 ) 2 ) 



i \2sis 2 m* - 2(si + a 2 )(si*l + ^2)"^ + si^si 2 - 2t 2 si + s 2 2 - 2s 2 ii + 2(ii 2 + t 2 2 ))^ 



+ m 2 (si - s 2 ) 2 + 2sis 2 (si + s 2 - t\ - t 2 ) - 2ss\s 2 



(B23) 



Here we have used the three body final state Mandelstam variables using the notation of II141L with the 
momenta assignments e~(p a ), e + (pb), f(pi), V(p 2 ), f(ps) and e is the electromagnetic coupling constant 
while Q is the charge of a particular quark. The corresponding production cross section at LEP is given by 



a(e+ e-^i'^Vff) 



1 



ds 1 ds 2 dtidt 2 \M\ 0[-Q] 



(B24) 



2 1( W J V-G 

where the physical region that determines the integration range satisfies the condition that the Gram deter- 
minant Q < 0. The Gram determinant for the 2 — > 3 production is explicitly given by 

Q = — \mls 2 - 2m 2 v s{sti + st 2 + si(s 2 - h) - s 2 t 2 + 2*it 2 ) (B25) 
lb 

+2s 1 (s(s 2 (t 1 + t 2 ) + h(t 2 - h)) + s 2 t 2 {h - 8 2 )) + (s(ti - t 2 ) + s 2 t 2 ) 2 + s?(s 2 - ti) 2 ] . 
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d. LHC gg — > SS production 

For a field in a flavor representation F, and color representation R, the cross section for g g — > S S is 

da dim(F) it a 2 



dt " 8 2 s 2 
Here the functions Fs are 



[d F 1 + C 2 F 2 + C 3 F 3 + C 4 F 4 ] . (B26) 



s — 2 m 2 s — 2m 



2 



F x = 2 + ^ + 



(«-t) 



2 



^2 = ^^ + 



9 1 9 ' 

j u — mi. 
(m 2 s - u)(Am 2 s - s) + (t - u)(s - 2m 2 ) /2 



2s : 



F 4 



s(t — m z s ) 



1 1 



(B27) 



(t — m 2 ) 2 (u — m 2 ) 2 J 
2(s - 2m 2 s ) 2 



(B28) 



(u — m 2 )(t — m 2 ) 
The color factors, for the R = (3, 6, 8) representations (G = 8) are 

d = C(R)dim(G)(2C 2 (R) - C 2 (G)/2) = (14/3,310/3, 108), 
C 2 = C(R)C 2 (G)dim{G) = (12,60,72), 
C 3 = [C 2 ( J R)] 2 dim( J R) = (16/3,200/3,72), 
C 4 = C(R)dim(G)(C 2 (R) - C 2 (G)/2) = (-2/3, 110/3, 36). 

Appendix C: Low Energy Constraints 

The main motivation for studying flavor symmetric extensions of the SM is to have an effective way of 
suppressing FCNCs allowing new states of fairly low mass. In this appendix we give a brief account of the 
constraints on our models from low energy data. We focus on processes that can be very restrictive, K — K, 
B s ,d ~ B s4 and D-D mixing. 

The effective weak Hamiltonian describing the B s — B s mixing in the SM is 

H*£ = VWVRG^^(sL^b L ) 2 , (CI) 



A 2 



where 



Asm = 2?F / -— = = 4.39 TeV. (C2) 

G F mw\/ &o{x t ) 

This is larger than m\y by an order of magnitude because SM contribution arises at one loop level from box 
diagrams with t and W exchanges. The factor r/w encodes the QCD corrections to the matching at the hard 
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scale nw ~ mf, while t]rg is the correction due to the running from fiw to In NDR at NLO they are 
r\w = 0.970, rfRG = 0.860 (using as(m z ) = 0.1184, fj,yy = m\y, ^ = m& = 4.2 GeV and m t = 164 
GeV for which So(x t ) = 2.33). The SM weak hamiltonians for the B d — Bp, K — K and D — D mixing 
follow from ( |C1| ) with the obvious replacements in quark flavors. In addition K and D mixing amplitudes 
receive large nonlocal contributions in the SM. 

NP contributions can modify the effective weak Hamiltonian to T~L^ + AT-L^. To make the comparison 
with the experiment easier we normalize the NP contributions to the SM predictions. For B ad mixing we 
thus define 

S ' d (B s , d \-H^\B s , d ) ■ ^ 

Note that a St d measures the deviation of NP weak phase from the SM one (for B s mixing the SM weak 
phase is very small, 0(1°), while for B d mixing it is 2/3 ~ 45°). Since the (real part) of {K^U^^K ) 
receives potentially large long distance contributions from AS = 1 operators, we normalize the NP matrix 
elements to the measured mass difference 

Amx ArriK 2rrtK 
Similarly there are potentially large long distance effects in D — D mixing amplitude, so that we define in 
the same way 

hoe ^ s ^ = ■ <°°' A < P I 5 °>. (C5, 
Amp ArriD 2m,£, 

Note that uk and <td measure the total weak phase of the NP contributions to the mixing, not just the 
deviation from the SM one. No NP in the mixing corresponds to h s = = kx = hp = 0. 

From the agreement of measured and predicted values of ex we obtain the bound at 95% C.L. 

\h K sin(2o-*r)| ^ 1-3 x 10~ 3 . (C6) 

Conservatively we use the the prediction \e K \ = (2.01^;^) • 10~ 3 from fl22l . which is in good agreement 
with the measured value |e^| = (2.229 ± 0.010) • 10 -3 (for a different treatment of lattice inputs see 
H140II ). The bound on hx cos 2a k is more uncertain because of potentially large long distance contributions. 
Assuming conservatively that Amu is saturated by NP gives \hx cos2ax\ 0.5. The CP violation in 
D ~ D mixing is well constrained and so \ho sin(2<7D)| ^5 0.3, at 95 % C.L. II 1421 (the NP contribution 
to 4>i2 is 2a d, while the SM CP violating contribution is negligible). Conservatively assuming that NP 
saturates the mass difference Amo we obtain fi£, < 0.5. 
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The NP contributions to the effective weak Hamiltonian in (CI I, AH^", are described by local operators 

A^Sf = Yl + ^Qi (C7) 

i i 

where the sum runs over the full set of dimension 6 local operators. For B s mixing these are I! 14311 

Ql = (sLj t Mb L )(si J 'y tl b L ), Q 2 ,4 = (sRb L )(s R>L b L ,R), Q3,5 = (»H,l&L,fl)> ( C8 ) 

Above Pl,r = (1 ± 7s)/2, while color indices are not displayed, if the summation is within brackets. The 



operators Q, L are obtained from Eq. (C8 1 by making a Pr •<-)• Pl replacement. Note that Q\ is the operator 



in the SM effective weak Hamiltonian Eq. ( |C1[ ). 



The SM and NP weak hamiltonians for the Bj, K and D mixing follow from \C\) and ( |C8| ) with the 
obvious replacements in quark flavors. In addition K and D mixing amplitudes receive large nonlocal 
contributions in the SM. We evaluate the matrix elements for the operators Qi in the same way as described 
in Eqs. (8)-(10) of M143II , including the numerical values, except for the updated values of fs d = 205 ± 12 
MeV, f Ba = 250 ± 12 MeV JT441 and the quark masses, for which we use the PDG 2010 values /3 = 
(21.8 ± 0.9)° and 7 = (67 ± 4)° ifUBI . The matrix elements for Q { are the same as for Qi since QCD 
interactions are parity conserving. 

1. Models I-IX 

The models I-IX give contributions to the effective Hamiltonian for B s mixing that at scale \x = raj take 
the form 

An ^' Bs = ^(VtVtbVZfv'VRG \b L -f8 L h^s L ] , (C9) 

or of its parity image: replacing 1)r^sr bR^^SR for SlT^sl &l7/xSl- Parity symmetry of the strong 
interactions implies that these operators have identical matrix elements between one particle states. The 



value of model dependent constant n s is listed in Table VIII while My stands for the average mass of the 
MFV-vector exchanged. The factor 7/ is a QCD correction that arises from running from the scale of the 
new physics, My, to the electroweak scale, My/- In calculating bounds we set 7/ = 1 incurring an error of 
order [a 5 (M w )/4vr] ln(M v /M w ) « 0.01 \n(M v /M w ). 

The contributions to B^, D and K mixing we define analogously (for simplicity of notation setting 
V'VRG = 1)> 

=^yt V M) 2 (d L ^b L ) 2 , AH^> K = ^(ytV ts V: d ?{d Ll »s L )\ (C10) 

A ^ NffP,D =^yl v ^\c L ^uLf. (en) 
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Case 


AC 

K s,d/C s ^ Kx/C Si0 Kd/C s _ 


Is,o 


m d R tyY D &uYl)d R + [ m d R (yY D AuA D Yl)d R + h.c] 

vldvUm +myl) 2 vlvlim + mvl) 2 o 


Hs,o 


m ur^YuAdY^ur + [772 urQ/YuAuAdY^ur + h.c] 

ylvlim + myl? 


IH S ,o 


rnQ L <y&u)Q L + v'iQl(YA d )Q l + [ V2 Ql(Y A V A D )Q L + h.c] 

{m + mvl) 2 vi vi 2 


IV S , 


VoQlVQl + ViQl^A v )Q l + ^Ql^A^Ql + [ V2 Q L ^ A V A D )Q L + h.c] 
(m + mvl? + A C c r)\ + A cc V? + Acc 


V s , 


Thd R Vd R + [md R Y{Y D AuYl)d R + h.c] + m d R (Y D A v Yl)Y (Y D A v Yl)d R 

,,2 „2 [ Vl 1 (»72-2r)J) 2 r]±7]3 ] 2 „ .2 \ {m-^lf , im+vl) 2 , V1V3] n 
Vs,dyb |_4r 2 + 12r 2 1 r 2 _ 4 J 4^24 [ 4r| 1 12r 2 + rf ] U 


VI S)0 


rjiURfuR + [tjzUrY (YuA D Y^)uR + h.c] + m u R (Y u AnY^)Y(Y u A D Y^)u R 

n n ii 2 n 2 | (m+vZ) 2 I 7?U73 

U U VuVc [ 4r§ 1 12rf + r 2 J 


VII S , 


?7i <iflf u R + 772 ^(YdAuY^)^ + h - c - 
Acc Acc Acc 


VIII S , 


ViQlYQl + [v2QlAuYQl + t)' 2 Q l AuYQl + h.c] + ^QlAuYAuQl + v' 3 QlA d YA d Q l 

V2 2 1 (')2- 2 ')2) 2 | VlV3 (Vi-mf , (m+nl) 2 | J?ll)3 (V 2 -n'2*) 2 1 (')2+'fe*) 2 1 ^l^S 

Irf 12r 2 1 r 2 _ 4 4r 2 1 12r 2 1 r 2 4r| 1 12r| ' rf~ 


IX Si0 


ViQlYQl + [v2QlAuYQl + V^QlAdYQl + h.c] + ■■■ 

vf 1 (%-2>j2) 2 , A ('J2-'/2) 2 | (»72+%) 2 , A W2-V2*) 2 I (%+')2*) 2 I A 

4>J + 12r| 1 4r 2 1 12r 2 1 fl CC 4r | 1 12r | + A CC 



TABLE VIII: The part of interaction Lagrangian AC relevant for meson mixing amplitudes and coefficients Hi in 
( |C9| l-( |ClT] l f° r eacn of the models I-IX. We use the abbreviations Au = YyYjj, Ad = Y d Yjj and r.; — my i / my for 
the ratios of vector meson masses in the flavor multiplet to the average vector meson mass my, while Acc stands for a 
1-loop contribution, a box diagram, from exchange of charged MFV-vectors; see text. The factors of C s — 1(C = |) 
must be included for color singlet (octet) vectors. Due to lack of space we display for model IX results only for leading 
complex parameters (the undisplayed leading real terms are the same as for model VIII). 



The coefficients n s ^d,K,D are collected in Table VIII In models V, VIII and IX several vector fields mediate 
mixing and we indicate their separate contributions. This requires we distinguish among their masses , with 
My then denoting their average mass in the defining formulae ( C9 >-( C 1 1 > for Kj. 
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Numerically, for Bd )S mixing contributions we have 



1 AIm _,„„ i^TeVV 2 



,4 

The contribution to the K mixing is 

"- i2 ™ " WOW^ig^ - al x {W)' ■ <Cl3) 

and for the D mixing 

hnf- = WWi^^^ » (0.3 • 10-, x KD e«V . (C14, 

In the numerics above we use f K = 156 MeV, B K = 0.79 |[T43l . J3f = 0.865 |[T43l . the averages of 
CKM elements from CKMfitter ICHEP2010 update 111451 . The standard CKM unitarity triangle angles 
have the values /3 = (21.8 ± 0.9)° and 7 = (67.2 ± 3.9)°. We have also included the effect of running 
from electroweak scale to mn f or D — D mixing, however, this is numerically unimportant and results in a 
percent level shift in ho- 

The flavor constraints thus require that (for My = 1 TeV) the Ki coefficients for the particular model 
are below k& < 0.02, kk ^ 10~ 2 (unless the phase of kk is finely tuned not to give a contribution to ek), 
and kd ~ 20. In order to explain the hint for nonzero phase of B s mixing, on the other hand, requires 
h s ~ 0.02 or h s ~ 0.1 for the two solutions, with appropriate weak phases. 

We make the following observations regarding models I-IX: 

• We can group the models into two categories. The first category form "universal" models, where 
the contributions to the mixing amplitudes are due to class- 1 operators in the notation of Ifl3l . For 
universal models we have 

K s = K d ~ K K ~ KD, (C15) 

where the last approximate equalities are valid, if all the couplings r/j are of the same size. In addition 
in y b limit k s4 = k k . 

• The second category form the "Yukawa suppressed" models, where the contributions to the mixing 
come from operators with additional Yukawa suppression (the so-called class-2 operators in the nota- 
tion of [ 13]). These are the models II S , > VI S)OJ that contribute only to D — D mixing and the models 

I s ,o, V S)0 , for which 

K d : k s : k k = y\ : y 2 s ■ (C16) 

yb 

(the relation to kk is only approximate for model V). 
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• Universal models can require small values for r/j ~ 0.1 • (My /I TeV) in order to avoid K — K 
and Bd — Bd mixing constraints. This is an order of magnitude smaller then the flavor-conserving 
couplings typically required to obtain large A f p B . This by itself is not a problem as the flavor violating 
couplings could be loop suppressed. 

• Yukawa suppressed models have no problems satisfying FCNC bounds. The typical sizes of contri- 
butions to hi are 

^ (2 . 10 - w g(^) 2 (^) 2 , 

where we have normalized y s ^ to values typical for large tan (3 (yb ~ 1). The contributions to B&, 
K and D mixing are completely negligible for these models. 

• Yukawa suppressed models I S)0 , V s>0 can potentially explain B s mixing anomaly, if yb ~ 0(1) and 
f]i ~ 0(1), while vectors need to have electroweak scale masses of no more than a few 100 GeV. 
For this a flavor diagonal weak phase is needed [92]. A flavor diagonal weak phase arises if the 
coefficients r\i are complex. For models I- VI the coefficients can be real only for terms that are 
higher order in Yukawas and therefore non-hermitian: 7/0,1,3 and r( x are all real, while r/2 can be 
complex. In models VII S all r/j can be complex. In models V, VIII and IX the contributions to the 
mixing carry a weak phase only if the vector mesons are not degenerate. 

Finally, in models IV, VII and IX we include a contribution Acc from charged vectors that produce 
mixing through a 1-loop box diagram, much like in the standard model. For example, for model IV S we 
find for the contribution to B s , Bd and K mixing 



128tt 2 yl 
where x tv = (m t /M v ) 2 and [146] 



1 "fon, \ 



&cc = r—-^ J iS (x tv ), (C21) 



S (x) 



19 1 3 1 

+ 



4 41 -x 2(l-x) 2 
For model VII S the B s mixing contribution is 



Q / \ 3 

6 I X 



2 



ln(x) . (C22) 



Acc = ^(y b y s ) 2 S (x tv ) [\ mm \ 2 + Re(^r/ 2 ) 2 ] , (C23) 
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0.2 



0.1 



^ My/in, 

FIG. 14: Correction /i, in units of ijq, to BB mixing in Model IV S from double exchange of charged currents as a 
function of vector mass in units of top mass. 

while the Bd and K mixing contributions are obtained through the replacements y s — > yd and y^ — > yd, 
respectively. Above, we have assumed degenerate MFV-vectors for simplicity. For IV Q and VII insert 




an additional factor of 11/6 in these formulae. Fig. 14 shows the correction h = h s = hd in model IV, 



resulting only from Acc> as a function of the top to charged vector mass ratio mt/My. The figure is 
for 770 = 1, but can easily be rescaled as indicated. The bound h < 0.2 gives My > 200 GeV (using 
mt = 164 GeV) for t]q = 1, increasing to My > 420 GeV for 7/0 = 2. In estimating this effect we have 
assumed the box diagram for these vectors is precisely as for the case of W's in the SM, because the GIM 
mechanism renders the diagram finite even in unitary gauge. 

In the D — D case, in contrast to the B — B case, the long distance contribution is believed to dominate 
the double charge current exchange amplitude. Since this is not computable, we settle for estimating Acc 
by comparing the perturbative box diagram to the corresponding one in the SM. Since all internal quarks 
are light the computation simplifies: for example, in model IV S the ratio is 

NPbox = (mV (MwV „ ( 19077Q V 

SM Box W \M v ) ~\My(GeV)J ' K } 

so that, for example, the new physics is a factor of 10 suppressed relative the SM for My/770 ^ 340 GeV. 

2. Models X and XI 

The vector models give a contribution to the effective Hamiltonian for B s mixing that at the scale /j ~ 
my takes the form 

An ^ ,Ba = w {y2tVtbK)2 (Q4 - Qs) ' (C25) 

where the operators (^4,5 are given in (|C8|). The contributions to Bd, K and D mixing are defined in the 



same way, with appropriate replacement of the CKM factors and quark fields, as in (C10 1, (Cll ). The above 
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Case 



X 



3.0 



AC 



K K 



Kd 



VidnVQl + mdnY^jQl + m d R (y D AuYl)yQi + m d R (Y D AuY^ )f A* V Q C L + h.c. 



Vs,dVb 



mvi i Viva 



l"i3 



V2V3 

r 3 i 



VsVd 



^ 



ViVi 



thurYQI + r) 2 u R yA* D Q c L + V3 u R (YuA D Y^)YQ c L + r, 4 u R (YuA D Y^)YA* D Ql + h.c. 



Xlq 



TABLE IX: The same as Table VIII but for models X and XI. We further defined rji = ?yi + 772 J/f + VsUtDb + ViUtVb' 



V2 = V2 + ViUt' V3 = 7 l3 + ViVt > an ^ Vi — Vi + ^22/6 + VsviVt + VAVbVt- The results are written in the limit of 
degenerate vectors mediating transitions between the first two generations. Above the index i stands for i = 1,2. 

expression is for color triplets, Xg, Xlg. The result for color sextets, Xq, Xlg, is obtained with a replacement 
(Q4 — Q5) — > — (Q4 + Q§)/2. The factors k« are collected in Table IX 



To obtain the low energy predictions, we include the RG running using the equations and numerical 
values for the matrix elements given in 111431 . As before, we use the updated values of CKM elements from 
CKMFitter, ICHEP2010 update HH. We obtain 



h,e 



i2<j a 



1.7(-2.1) x y? 



4 k*s ( y. 



Vs V0.02 



1 TeV 

my 



h d e 



i2a d 



0.09(-0.11) x y? 



A K d ( Vd 



Vd 



IO- 3 



i2a K 



5(-4) • 10" J x e up y 4 t 



i2/3 4 K K ( jy_ 



VdVs 



lTeVV 
my J 

Vd 



0.02/ VIO- 3 

Vu 



ITeV 



my 



Vc 



h D e<*"> * 1.4(-1.5) • 10- 8 x e^yt^ y , , r 

y u y c V10 V \5 • 10 6 



ITeV 

my 



(C26) 
(C27) 
(C28) 
(C29) 



where the first values are for color triplet models and in the brackets for the color sextet models. 
We make the following observations: 

• Models X and XI are of the "Yukawa suppressed" category, but the class-2 operators that they give 
rise to are linear in y s , y d suppression, which is quite atypical. Commonly the suppression in class-2 
operators is quadratic in light Yukawas. The ratios of Ki for models X are 



K d : k s : k k = y d : y s 
where the ratio with kk is only approximate. 



VsVd 
Vb 



(C30) 
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• Models Xg 6 offer the possibility to explain A l p B and B s mixing anomaly at the same time, if fla- 
vor violating couplings rji are a factor few smaller than the corresponding couplings in the flavor 
conserving part of the Lagrangian. 



3. Models Si_iv 



The models Si_rv give contributions to B s mixing of the form 



*H™' B ' = ~ 2 (y 2 tV tb V: s fQ„ (C31) 



where the operators Q4 5 are given in (C8 1. The contributions to B d , K and D mixing are defined in the 



same way, with appropriate replacement of the CKM factors and quark fields, as in (CIO I, ( Cll[ ). The above 



expression is for models Si and Sin, while for models Sn and Srv the result is obtained with a replacement 
Qi — > —Q4/6 + Q5/2. The factors Ki are collected in Table |X) Numerically these models are very similar 
to the MFV models X and XI, 

-h a e»* * 1.5(0.05) x yf^ (J£-) (C32) 

y s V0.02/ \ ms J 

-had*" * 7(0.3) • 10- 2 x yf^ ( W ) (1™) 2 , (C33) 

-h K e">« - 3.6(-0.1) • 10- 3 x e* 2 V— (^) (£*) ( — Y , (C34) 

y d y s V0.02/ V10 A ) \ m s ) 

-hoe*"* * 1.1(0.01) • 10- 8 x e^y*— ( ( — ^-V) ( ^^-) ' , (C35) 
V ' b y u y c V10" 5 A5-10- 3 / V rn s J 

where the values in the brackets are for the models Sn and Srv, f° r which numerical cancellation between 
the matrix elements of the two operators occurs. 



4. Models S v xiv 

Of the remaining models only models Svi,viii,xiv lead to tree level FCNCs. The models Sv,vii,ix,xi 
and Sxin do not lead to tree level FCNCs because of its color triplet structure (and thus antisymmetric 
contraction of indices). The models Sx and Sxn lead to charge currents and thus give FCNCs at loop level 
only. 

We next focus on the models that lead to FCNCs at tree level. They give contributions to B s mixing of 
the form 

^C^-SfeWQi, (C36) 
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Case 






AC 

K K 






rj u R SQ L -+ 


■tjxUrSAdQl - 


hV2U R (YuA D Y^)SQ L 


+ muR(YuA D Y^)SA D Q L + h.c. 


i, ii 




















tf„yc[2^2^2l2al + JiSgi!] 


III, IV 


r] d R SQ L -f 


- ijxcIrSAuQl - 
_l_ V3Vp _|_ mnl 1 

r 3i r 33 J 




+ VsdRiYoAuY^SAuQL + h.c. 
\mmil] o 



TABLE X: The same as Table VIII but for models Si_iv- We also use the abbreviations ?/o = Vo + ViVt > ?7i = 
'/i + ihUbUt ' ^2 = ? /2 + '732/?- The results are written in the limit of degenerate scalars mediating transitions between 
the first two generations. Above the index i stands for i = 1, 2. In all the models there is also a 1-loop contribution 
Ace ( not displayed). 

where the operator Q\ is replaced by Q\ for model Sxiv- The contributions to K and D mixing are 



defined in the same way, with appropriate replacement of the CKM factors and quark fields, as in (CIO), 



(JCTTJ). The factors k- l are collected in Table |XI| Numerically the analysis is exactly the same as for MFV 
vector models Si_x- The Yukawa suppressed (class-2 operator) models are Syi and Sym, while Sxiv 1S an 
example of the universal model (class- 1 operator model). 



5. Models S H ,8 
The models Sh,8 give contributions to B s mixing of the form 

AH^' Ba =k s (— 2 \- ) {y 2 t V th V t * s ) 2 {C 2 Q2 + C 3 Q 3 ), (C37) 

l m s, m sJ 



where the operators Q 2t3 are given in ( |C8| ), ms 1 2 are the masses of CP even and odd neutral scalars, and 



C 2 = 1/2(1/4), C 3 = 0(-l/12) for models S H (S 8 ). The contributions to B d , K and D mixing are 



defined in the same way, with appropriate replacement of the CKM factors and quark fields, as in ( |C10[ ), 
( |C11| ). Numerically (taking ms 2 S> ms 1 limit for simplicity), 

ITeVV 

m s 



h de ^s. d „ _ 15(5) x 4 (C38) 

"Si / 



h K ^« - l.l(-0.3) • 10- 2 x fJ^Y (1™L) , ( C 39) 

yj V0.02/ V m Sl ) 

h D fP"° ~ -1.1(0.3) • 10- 6 x e^yt^ (^^Y (—)\ (C40) 

y£ V5 • 10 6 J V m Si / 
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Case 


K s,tZ 






VI 


t]qu r du r -+- rjiu R D(iu ^D-f u) u R 




r l2U R {Iu'-^D^u) ^(^U^D^jjfUR -t- n.c. 










2n( mV [* + VoV ' 2 *)y 2 u y 2 c 


VIII 


jjo^dfl + rjid R S(Yi)A.i/Y^ ) )dji H 


- V ' 1 d R (Y D A u Yl) T Sd R + 


Tl20. R \IDi-^U Y jj) *J\i Dl-±U * D) a R 1 n - c - 






y 2 d y 2 s 2$t(mrii* + V0V2 


) 


XIV 




-^QlAjjSAuQL + h.c. 




r23?r;ir;j* J7 J)J . i)2l *l 




2SR(^77i* + r ?0 r?|) 



TABLE XI: The same as Table VIII but for models Sy-xiv- Only models that lead to tree level FCNCs are shown. 
The results are written in the limit of degenerate scalars mediating transitions between the first two generations. Above 
the index i stands for i = 1, 2. For model Sxiv the terms that contribute to D — D mixing have not been written out 
explicitly in AC. They can be obtained with Au —> Ad and have coupling fji. 



Case 



AC 

K K 



H,8 



ilodRY D SQ L + tj'qUrYuSQl + riidR(Y D Au)SQ L + rj[u R (Y u A D )SQ L + h.c 

yWi + Acc yh? + A 



vWi + A cc 



cc 



TABLE XII: The same as Table VIII but for models Sh.s- The 1-loop contributions Acc arise from charge currents. 



where the values in the brackets are for model Ss- 
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